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EQUATION WITHOUT ANGULAR CUT-OFF 
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Abstract. We prove the existence and exponential decay of global in time 
strong solutions to the Boltzmann equation without any angular cut-off, i.e., 
for long-range interactions. We consider perturbations of the Maxwellian equi- 
librium states and include the physical cross-sections arising from an inverse- 
power intermolecular potential r _ ' p_1 ' with p > 3, and more generally, the 
full range of angular singularities s = u/2 € (0,1). These appear to be the 
first unique global solutions to this fundamentally important model, which 
grants a basic example where a range of geometric fractional derivatives oc- 
cur in a physical model of the natural world. Our methods provide a new 
understanding of the effects of grazing collisions in the Boltzmann theory. 
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1. Introduction, main theorem, and historical remarks 

In 1872, Boltzmann was able to derive an equation which accurately models the 
dynamics of a dilute gas; it has since become a cornerstone of statistical physics [15, 
16,22,27,28,51]. There are many useful mathematical theories of global solutions 
for the Boltzmann equation, and we will start off by mentioning a brief few. In 
1933, Carleman [14] proved existence and uniqueness of the spatially homogeneous 
problem with radial initial data. For the spatially dependent theories, it was Ukai 
[48] in 1974 who proved the existence of global classical solutions with close to 
equilibrium initial data. Ten years later, Illner-Shinbrot [34] found unique global 
mild solutions with near vacuum data. Then in 1989, the work of DiPerna-Lions 
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[22] established global renormalized weak solutions for initial data without a size 
restriction. We also mention recent methods introduced in the linearized regime 
by Guo [32] in 2003 and Liu- Yang- Yu [38] in 2004. All of these methods and their 
generalizations apply to hard sphere particles or soft particle interactions in which 
there is a non-physical cut-off of an inherently nonintegrable angular singularity. 

When the physically relevant effects of these angular singularities are not cut-off, 
the only global spatially dependent theory we are aware of is the remarkable paper 
by Alexandre- Villani [5] from 2002, which proves the existence of DiPerna-Lions 
renormalized weak solutions [22] if one can add to the equation a non-negative 
defect measure. It is illustrated therein that the mass conservation they prove 
would imply this defect measure was zero if the solutions were sufficiently regular. 
At the moment this defect measure appears difficult to characterize [5, Appendix]. 

Despite the well-known physical and mathematical importance of this problem, 
it is perhaps the last remaining physically relevant case in the Boltzmann theory in 
which, as far as we know, there is no theory of global in time strong solutions for 
spatially dependent initial data of any kind. This issue is mentioned particularly 
in Villani [51]. The results herein prove the existence of such solutions for cross 
sections arising from an an inverse-power intermolecular potential r~( p ~ is> with 
p > 3 and, more generally, for the full range of angular singularities. 

Let us now give a detailed explanation. We study the Boltzmann equation 

(1) ^+v-V x F = Q(F,F), 

where the unknown F(t,x,v) is a nonnegative function. For each time t > 0, 
F(t, •, ■) represents the density function of particles in the phase space; some may 
call F the empirical measure. The spatial coordinates we consider are x € T 3 , and 
the velocities are v € K 3 . The Boltzmann collision operator Q is a bilinear operator 
which acts only on the velocity variables F(v) and is local in (i, x), as 

Q(F,F)(v)= I dv, f da B(\v-v*\,a)[FlF' -F*F]. 

JR 3 JS 2 

Here we are using the standard shorthand F = F(v), F* = F' — F(v'), 

Fl = In this expression, v' , v'^ and v, are the velocities of a pair of particles 

before and after collision, and are connected through the following formulas 

, v + v* . \v-v*\ , v + v* \v-v*\ ^ S 2 

There are other ways to represent Q which result from alternate choices for the 
parameterization of the set of solutions to the physical law of elastic collisions: 

(6) M 2 + M 2 = M 2 + K| 2 . 

We specifically discuss Carleman-type representations in the appendix of this paper. 

The Boltzmann collision kernel B(\v — v*\,cr) for a monatomic gas is, on physical 
grounds, a non-negative function which only depends on the relative velocity \v— u*| 
and on the deviation angle 9 through cos# = (k,a) where k = (v — v*)/\v — v*| and 
(■, •) is the usual scalar product in R 3 . Without loss of generality we may assume 
that B{\v — v*\,o-) is supported on (k,a) > 0, i.e. < 9 < ^. Otherwise we can 
reduce to this situation with the following "symmetrization" : 

B(\v-v*\,(t) = [B(\v-v*\,<T)+B(\v-v*\,-a)]l( kt>T )> . 
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Above and generally, 1a is the usual indicator function of the set A. 

1.1. The Collision Kernel. Our assumptions are as follows: 

• We suppose that B(\v — v*\, a) takes product form in its arguments as 

B(\v-v*\,a) = $(|u-t>*|)6(cos0). 

In general both b and <!> are non-negative functions. 

• The angular function t h-> b(t) is non-locally integrable; for Cb > it satisfies 

(4) _^< sin 06(cos0)<_I_ ,€(0,1), V0e(o,| 

• The kinetic factor z H> 3>(|z|) satisfies for some C$ > 

(5) $(> - v*\) = C$|w-w*| 7 , 7 > -min{2s,3/2}. 
Notice that 7 + 2s > ensures a spectral gap from [42] . 

Our main physical motivation is derived from of particles interacting according 
to a spherical intcrmolecular repulsive potential of the form 

0( r ) =r -O-i) ! pe (2,+oo). 

For these potentials, Maxwell in 1867 showed that the kernel B can be computed. 
It satisfies the conditions above with 7 = (p — 5)/(p— 1) and s = l/(p— 1); see for 
instance [15,16,51]. In this situation 

p-3 

7 + 2s = > 0, or p>3. 

p- 1 

Thus all of the conditions in (4) and (5) hold with — 1 < 7 < 1 and < s < 1/2. 
Some authors use the notation v/2 = s € (0, 1) which is equivalent to our own. 

Many research results on the non cut-off Boltzmann equation consider regular- 
ized kinetic factors, which means that $(\v — u*|) = C$(w — v*) 7 . In this situation 
our results apply easily to any s € (0, 1) and any 7 > —2s; we can also handle these 
exponents in the physical kinetic factor from (5). We do elect not to record these 
details herein since it would add unnecessary technical complexity that does not 
directly relate to the main goal of illustrating our new methods. We can further 
prove our main result for the full physical range: 7 + 2s > —1 and p > 2; this 
will be the content of a forthcoming work [29]. However, in the context of our 
main theorem and [29], the essential mathematical difficulties associated with the 
angular singularities (4) are resolved in the present result. 

We will study the linearization of (1) around the Maxwellian equilibrium state 



(6) F(t,x,v) = fi(v) + y/fi(v)f(t,x,v), 
where without loss of generality 

fi(v) = (2^)- 3 / 2 e-^ 2 / 2 . 

We will also suppose without restriction that the mass, momentum, and energy 
conservation laws for the perturbation f(t,x,v) hold for all t > as 

(7) [ dx dv ( v ] y/n{v) f{t,x,v) = 0. 

^T|XB3 \M 2 / 

This condition should be satisfied initially, and then will continue to be satisfied 
for a suitably strong solution. Our main interest is in global strong solutions to 
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the Boltzmann equation (1) which are perturbations of the Maxwellian equilibrium 
states (6) for the long-range collision kernels (5) and (4). 

As will be seen, our solution to this problem rests heavily on our introduction of 
the following new weighted geometric fractional Sobolev space: 

N^^{feL 2 (Rl) : |JV„<cx>}, 
where for v G R 3 we define v= (u, ^ |i>| 2 ) G R 4 , and then specify the norm by 

(8) +/ * / dv>((v)(v'))^ ^"g liWKi. 

This space includes the weighted L 2 space given by 

l/l 2 L2 d = f / dv ( V y +2s \f(v)\ 2 . 

Note that the weight is as usual: (v) = yT+ |i>| 2 . The inclusion of the quadratic 
difference |i>| 2 — \v'\ 2 in the fractional kernel will be of great importance; it is not 
a lower order term. The rest of our notation is defined below. 

The sharp space iV s ' 7 is equivalent to a weighted, nonisotropic Sobolev norm, a 
feature which was conjectured in [42] . Precisely, if R 3 is identified with a paraboloid 
in R 4 by means of the mapping v i-> (v, ^\v\ 2 ) and Ap is defined to be the Laplacian 
on the paraboloid induced by the Euclidean metric on R 4 then 



mk, « / dv(vy +2s \(i-A P )if(v) 
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We will however omit the proof of this characterization as it has no direct role in 
establishing our results. With this, we may now to state our main result as follows: 

Theorem 1. (Main Theorem). Fix N > 3, the total number of spatial derivatives. 
Choose fo = fo(x,v) G L 2 (R 3 ; iJ^CF 3 )) in (6) which satisfies (7) initially. There 
is an r)o > such that if ||/o||_L 2 (R3;ir N (T 3 )) — VO; then there exists a unique global 
strong solution to the Boltzmann equation (1), in the form (6), which satisfies 

f(t,x,v) G Lr([0,oo);iX w (T 3 x R 3 ,)) n i?([0, oo); N^H^ x R 3 )). 
Moreover, we have exponential decay to equilibrium. For some fixed A > 0, 

\\f\\L 2 (Rl;H N (Tl))(t) ^ e _A * H/oll Z, 2 (Kg;/f W (T|)) - 

We also have positivity, i.e. F = fi + WJIf > if Fq = fj, + i/fifa > 0. 

Grad proposed [28] in 1963 the angular cut-off which requires that &(cos 9) be 
bounded. Grad also pointed out that many cut-offs are possible. In particular, the 
following less stringent L X (S 2 ) cut-off has become fashionable 1 

da b((k, a)) < oo. 

s 2 

These types of truncations have been widely accepted, and have now influenced 
several decades of mathematical progress on the Boltzmann equation. For the 
intermolecular repulsive potentials previously discussed, the cut-off theory only 
applies physically in the limit when p — > oo, which represents Hard-Sphere particles. 

These cut-off assumptions were originally believed to not change the essential 
nature of solutions to the equation. It has been argued by physicists, see [51], that 



1 Note that this L 1 (§ 2 ) cut-off was already implicitly used in 1954 by Morgenstern [39]. 
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the important properties of the Boltzmann equation are not particularly sensitive 
to the dependence of the collision kernel upon the deviation angle, 9. Our theorem 
above shows at the mcsoscopic level of Boltzmann that this physical heuristic misses 
a strong dependence of the solutions on the angular singularity; specifically, there 
is a gain of velocity regularity and velocity moments globally in time. 

Results regarding these types of smoothing effects can be seen in the pioneering 
work of Dcsvillettcs [18] from 1995, which applied to simplified models such as 
the Kac equation. In a very recent preprint from 2009, perhaps the first local 
existence theorem for large data and moderate angular singularities < s < 1/2 
has been shown by Alexandre-Morimoto-Ukai-Xu-Yang [8] for initial data which is 
somewhat smoother than ours, they also show the C^° x v regularizing effect. By 
contrast, under the angular cut-off assumption the solution is known to have the 
same regularity in a Sobolev space as the initial data. These results go back to 
Boudin-Dcsvillcttcs [13] in 2000 for solutions near Vacuum, and this same effect 
has been recently shown in the near Maxwcllian regime [12,24]. 



As a result of the fact that the angular singularity (4) is not integrable on a 
sphere, it has been conjectured in numerous works that the nonlinear collision 
operator should behave like a fractional (flat) Laplacian in the velocity variable v: 



Our precise work at the linearized level shows that this conjecture is not the whole 
story. Certainly we see that there is a smoothing effect globally in time. However 
for our results the most useful intuitive point of view is to think of the collision 
operator as a fractional Laplacian on a manifold, and this manifold depends in an 
essential way on the collisional geometry. 

By comparison, the Landau equation, derived in 1936, is probably the closest 
analog we have to the Boltzmann collision operator for long-range interactions; 
however the Landau operator involves regular partial derivatives rather than frac- 
tional derivatives and for that reason may be somewhat more understandable at 
first. Landau's equation is obtained as the limiting system when p — > 2 in the 
inverse power law potential, the collision operator can be shown to satisfy [51]: 



Notice there is a metric of sorts in this case-in the -which depends in an essential 
way on your unknown solution F. Even in the simplest case when your unknown is 
the steady state, F = n(v), this oy weights more heavily angular derivatives [30]. 
In the general case a^- is known to be degenerate and not comparable to 5ij at 
infinity, see e.g. [17]. 

For this paper, the basic new understanding which enabled our progress was to 
identify that the fractional differentiation effects induced by the linearized Boltz- 
mann collision operator are taking place on a paraboloid in R 4 . While we do not 
directly identify dependence of geometry on the function F itself in our formula- 
tion, it seems only reasonable to suspect that such dependence may be relevant to 
future work. Before reviewing the details of our proof, we will mention past works. 

1.2. Review of the non cut-off theory. As above, it makes good sense to briefly 
review a few results for the Landau equation, which corresponds to the grazing 
collision limit s — > 1 of the Boltzmann equation for long-range interactions (see 



Q(F,F) 



Cp{— A V ) S F + lower order terms. 
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[6] and the references therein). For the spatially homogeneous case with hard 
potentials (roughly, replace l/|i>| above with |i;| 7+2 for 7 > 0), global existence of 
unique weak solutions and the instantaneous smoothing effect was shown for the 
first time by Desvillettes and Villani [21] for a large class of initial data in the year 
2000. Then Guo [30] in 2002 established the existence of classical solutions for 
the spatially dependent case with the physical Coulombian interactions (p = 2) for 
smooth near Maxwellian initial data in a periodic box. Guo's solutions were recently 
shown to experience instantaneous regularization in [17]. For further results in these 
directions we refer to the references in [17]. 

Due to length constraints, it is unfortunately not possible to give an exhaustive 
review of results in the non cut-off theory. We will however try to mention a sample 
of results. In the case of Maxwell molecules, e.g. p = 5, it is remarkable that the 
spectrum and cigenfunctions of the linearized collision kernel can be computed 
explicitly as was performed in a classical paper [53] from 1952. This was later 
simplified by Bobylev [11]; in this work the Fourier transform of the Boltzmann 
collision operator was shown to have an elegant form, which is now called Bobylev's 
identity and has found widespread utility. Pao [44] in 1974 used the early techniques 
of pseudodiffcrential operators and Bcsscl functions to study the spectral properties 
of the linearized operator for general inverse-power intermolecular potentials with 
p > 3. In 1981-82, Arkeryd in [9, 10] proved the existence of weak solutions to 
the spatially homogeneous Boltzmann equation when < s < h. Then Ukai [49] 
in 1984 obtained a local Cauchy-Kovalevskaya type theorem in a function space 
which was analytic in x and Gevrey in v; this work applied to the moderate angular 
singularities < s < ^, the convergence of Grad's cutoff approximation and the 
positivity of solutions was also established. In 1998 Villani [50] introduced the new 
spatially homogeneous weak H-Solution formulation which can handle all physically 
meaningful interactions with p > 2 and more generally; we refer the reader to this 
paper for a fairly exhaustive list of references up to 1998. 

In 1998 Lions proved a functional inequality [37] which bounds below the "en- 
tropy dissipation" by an isotropic Sobolev norm H" up to lower order terms, for a 
certain range of a. Then in the work of Alexandre-Dcsvillettcs- Villani- Wennberg [4] 
from 2000, this entropy dissipation smoothing estimate was obtained in the isotropic 
space H% with the optimal exponent s. This work further introduced elegant formu- 
las, such as the cancellation lemma and isotropic sub-elliptic coercivity estimates 
using the Fourier transform. This was in several ways the starting point of the 
modern theory of grazing collisions. Subsequent results of Alexandre- Villani [1,5] 
developed a rcnormalizcd DiPcrna-Lions theory of weak solutions with defect mea- 
sure, and established the appearance of strong compactness. Desvillettes- Wennberg 
[20] further demonstrated that solutions to the spatially homogeneous Boltzmann 
equation for regularized hard potentials enter the Schwartz space instantaneously. 
And recently in 2009 Desvillettes-Mouhot [19] proved the uniqueness of spatially 
homogeneous strong solutions for the full range of angular singularities s € (0, 1), 
and they have shown existence for moderate angular singularities s £ (0, 1/2). 

Broadly speaking, the approaches outlined above use the Fourier transform to in- 
terpret the fractional differentiation effects isotropically. Other interesting methods 
have been introduced to further study the fractional differentiation effects isotropi- 
cally, using more involved methods from pseudodiffcrential operators and harmonic 
analysis. In particular, some uncertainty principles in the framework of Feffcrman 
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[25] were introduced in 2008 by Alexandre-Morimoto-Ukai-Xu-Yang [7] to study 
smoothing effects of the Boltzmann equation with non-cutoff cross sections. Im- 
portant methods introduced in that paper [7], as well as Alexandre [3], Morimoto- 
Ukai-Xu-Yang [40] , and the references therein, establish the hypoellipticity of the 
Boltzmann operator. They develop methods for estimating the commutators be- 
tween the Boltzmann collision operator and some weighted pscudodiffercntial oper- 
ators. And they sharpen some of the isotropic coercivity and upper bound estimates 
for the Boltzmann collision operator. Recently these authors have new preprint [8] 
with a large data local existence theorem as previously mentioned, note also the 
earlier local existence theorem between two moving Maxwellians [2] of Alexandre. 

We also mention the linear isotropic coercivity estimates from [41,42]. In par- 
ticular Mouhot-Strain [42] identified the sharp weight, 7 + 2s, in our norm (8). 

During the course of the proof of our main Theorem 1, we develop a new point 
of view and a new set of tools for the long-range interactions, which we believe 
have implications for a variety of future results both in the pcrturbative regime 
and maybe even beyond it. We do not use any of the major non cut-off techniques 
described above, most of which are designed around the Fourier transform and 
isotropic estimates. Moreover, we do not study the Fourier transform of the collision 
operator at all. 

From the standpoint of harmonic analysis, the estimates we make for the bilinear 
operator (10) arising from our ansatz (6) fall well outside the scope of standard 
theorems. The operator and its associated trilinear form may be expressed in terms 
of Fourier transforms as a trilinear paraproduct; such objects have been the subject 
of recent work of Muscalu, Pipher, Tao, and Thiele [43] and are known to be very 
difficult to study in general. Known results for such objects fail to apply in our case 
because of the loss of derivatives (meaning that at least two of the three functions 
g, h, and / must belong to some Sobolev space with a positive degree of smoothness). 
Moreover, routine modifications of known results (for example, composing with 
fractional integration to compensate for the loss of derivatives) also fail because of 
the presence of a fundamentally non-Euclidean geometry, namely, the geometry on 
the paraboloid. This nontrivial geometry essentially renders any technique based 
on the Fourier transform difficult to use herein. Instead, we base our approach 
on the generalized Littlewood-Paley theory developed by Stein [45]. Rather than 
directly using semigroup theory, however, we opt for a more geometric approach, as 
was taken, for example, by Klainerman and Rodnianski [36]. Since the underlying 
geometry we identify is explicit, we are able to make substantial simplifications over 
both of these earlier works by restricting attention to the particular case of interest. 

1.3. Possibilities for the future, and extensions. We believe that our general 
methods and point of view can be useful in making further progress on multiple 
fronts in the non cut-off theory. Herein we list some of those. 

We will soon address the generalization to the very soft potentials and 7 + 2s < 
in a subsequent paper [29]. Furthermore, the generalization to the whole space case 
can be handled by combining our estimates with the existing cut-off technology 
in the whole space. It would also be interesting to prove the instantaneous C°° 
regularizing effect for our solutions in the spirit of [17]. 

We are hopeful that the estimates we prove can help to resolve the existence 
question for the Vlasov-Maxwell-Boltzmann system without angular cut-off; notice 
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at the moment the theory here is limited to the hard-sphere interactions [32]. It 
would be interesting to study the smoothing effect, which could potentially improve 
dramatically the strong assumptions in the current time decay theory [47]. 

For spatially homogeneous solutions, our results provide additional information 
for the high singularities s > 1/2, in which otherwise there is currently no existence 
theory for strong solutions [19,40]. The methods and point of view in this paper 
may help to treat the high singularities with large spatially homogeneous data. 

Lastly, it would be quite important to work with the estimates herein and in [29] 
to justify rigorously the validity of Landau approximation near Maxwellian. 

In the next Section 2, we linearize the Boltzmann equation (1) around the pertur- 
bation (6) and then explain the sharp space associated with the linearized collision 
operator. We further define all the relevant notation and formulate and discuss the 
main velocity fractional derivative estimates. Then we describe several key new 
ideas in our proof, and outline the rest of the article. 

2. Notation, reformulation, the main estimates, and our strategy 

Throughout this section we will define the relevant notation for the problem. We 
also reformulate the problem in terms of the equation (9) for the perturbation (6). 
Last and perhaps most importantly, we will state and explain our main estimates 
towards the end of this section. 

Throughout this paper, the notation A < B will mean that an a positive constant 
C exists such that A < CB holds uniformly over the range of parameters which 
are present in the inequality (and that the precise magnitude of the constant is 
irrelevant). In particular, whenever either A or B involves a function space norm, 
it will be implicit that the constant is uniform over all elements of the relevant space 
unless explicitly stated otherwise. The notation B > A is equivalent to A < B, and 
Ak B means that both A < B and B < A. 

2.1. Reformulation. We linearize the Boltzmann equation (1) around the pertur- 
bation (6). This grants an equation for the perturbation f(t,x,v) as 

(9) d t f + vW x f + L(f) = T(f,f), f(0,x,v) = f (x,v), 
where the linearized Boltzmann operator L is given by 

L(g) = a^ 1/2 Q(m, y/Hg) n~ 1/2 Q(VJl9, m) 

= / dv* da B(\v - u*|,cos0) [g*M + gM* - g'^M' - g'M'*} M*, 
Jr 3 Js 2 

and the bilinear operator T is given by 

(10) T(g,h) = M~ 1/2 Q(VW,\//^) = / dv* [ daBM^ti - g*h). 

JR 3 JR 3 

In both definitions, we take 

M(v) = ^Kv) = (27r)- 3 / 4 e~M 2 / 4 . 

When convenient, we will without loss of generality abuse notation and neglect the 
constant (27r)~ 3 / 4 in the definition of M. Finally, we note that 



(11) 



L(g)^~r(M,g)-T(g,M). 
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This reformulation shows that it is fundamentally important to obtain favorable 
estimates for the bilinear operator T. 

We expand the main term of the linearized Boltzmann operator as 



r(M, 5 ) 



/ dv* [ da B(\v - v*\,cos6) [M'J - M*g] M*. 
Jr 3 Js 2 

We will now split this in parts whilst preserving the cancellations as follows: 
T{M,g) = [ dv* [ da B [g' - g) M'.M. - v{v) g(v), 

JR 3 JS 2 

where 



v{v) = I dv* / da B(\v - u„|,cos0) (M* - M'*)M*. 

JR 3 JS 2 

The first piece above satisfies a favorable identity. This will follow from the pre-post 
collisional change of variables as 

-/ dv I dv* [ da B(\v-v*\,cos6)(g' -g)hMiM* 
Jr 3 Jr 3 Js 2 



f dv 1 


dv* 1 


/R 3 Jl 


i 3 Js 



-If dv f dv* [ da B {g - g')h! M*M'* 
^ Jr 3 Jr 3 Js 2 

= \ [ dv I dv* I da B (g' - g)(ti - h)M'*M*. 
2 Jr 3 Jr 3 Js 2 

This shows the crucial point that herein is contained a Hilbert space structure. For 
the weight, it has been shown that we can make the splitting 

v{v) = v(v) + v K {v), 

where under only (4) and (5) it can be seen that 

v{v) w (v) 7+2s , and \v K {v)\<{vy. 

These estimates were established by Pao [44, p. 568 eq. (65), (66)] when studying the 
eigenvalue problem. Pao reduced this question to the known asymptotic behavior of 
confluent hypergeometric functions, making use of the general addition theorems. 

We further decompose L = N + K. Here N is the "norm part" and K will be 
seen as the "compact part." The norm part is then written as 

(12) Ng = -V{M,f)-u K {v)f = -( dv* [ da B(g' - g)M'*M* + v{v)g{v). 

Jr 3 Js 2 

We will use (•,•) to denote the standard L 2 (M.^) inner product. Then, with the 
previous calculations, this norm piece satisfies the following identity: 



{Ng,g) = \! dv f dv* f daB(g' -g) 2 M'*M*+ f dv v(v) \g(v)\ 2 . 
1 Jm? jr 3 Js 2 Jr 3 



As a result, in the following we will use the anisotropic fractional semi- norm 

(13) \g\% = 



1 dv 1 


dv* 1 


Ir 3 J 


R 3 Js 
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For the second part of (Ng,g) we recall the norm 1/1^2 ^ defined below equation 
(8). These two quantities will define our designer norm, it is sharp for the linearized 
operator. We also record here the definition of the "compact piece" K: 

(14) Kg = v K (v)g - T(g, M) = v K (v)g - [ dv* [ daBM^M' - g*M). 

This is our main splitting of the linearized operator. 

To simplify many statements below, we will use the following (77, ^)-norm: 

1/1^=/ dvimfivivy^ + S- 1 ), r,,S>0. 

This norm perhaps requires some explanation. It will serve to unify several different 
desirable inequalities later. The first useful feature of these norms is to observe that 

( 15 ) i?f \g\s,ri\h\r,,s = \g\ L 2\h\ L 2 + \g\ L * ,,\h\v>- 
The second desirable inequality to be made is that 

\f\l n <m\l ' a dv\f(v)\ 2 , \/ v >o. 

This holds because 7 + 2s > so the weight 1 is ultimately bounded by any small 
constant times (u) 7+2;i provided \v\ is large enough. The right-hand side of this 
inequality is precisely the sort of norm needed for estimates of the "compact" part. 

All of the above L 2 (R 3 )-based norms are exclusively in the velocity variables. 
We will define analogous £ 2 (T 3 xlj) norms in both space and velocity variables 
by replacing | • | with || ■ | and keeping the same norm subscript. We also use (-, •) 
to be the usual £ 2 (T 3 x R 3 ) inner product. In particular 

= II |ft|jV«.7 || £ 2( T 3), \\h\\ 2 ^ S = || \h\ Vl g || £ 2( T |) ■ 

Now, for a multi- index a = (or, a 2 , a 3 ), we will use the spatial derivatives 

d a = d al d a2 d a3 

We will also take L 2 (R 3 V ;H N (T 3 )) with N > 3 spatial derivatives to be 

H' l lli 2 (R3;.H" JV (Ti)) = X/ H^ Q ^Hi 2 (R?xT|)- 
\a\<N 

For brevity, in each of these norms we sometimes use the notation L 2 , L 2 L 2 , L 2 H^ , 
etc, to denote these spaces without confusion. For example, we will sometimes write 
||/l|li2( R 3 ;J? jv( T 3)) = H^lli^N and use other norms such as ||/i|| 2 2 = |H| 2 2 (T 3)- 

2.2. Main Estimates. We will prove all of our estimates for functions in the 
Schwartz space, 6>(R 3 ), which is the well-known space of real valued C°°(M 3 ) func- 
tions all of whose derivatives decay at infinity faster than the reciprocal of any 
polynomial. Note that the Schwartz functions are dense in the nonisotropic space 
N sn , and the proof of this fact is easily reduced to the analogous one for Euclidean 
Sobolcv spaces by means of the partition of unity as constructed, for example, in 
section 6.2. Moreover, in all of our estimates, none of the constants that come 
up will depend on the regularity of the functions that we are estimating. Thus 
using routine density arguments, our estimates will apply to any function in iV s-7 
or whatever the appropriate function space happens to be for a particular estimate. 
Our first non-linear estimate is the following: 
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Lemma 2. For the non-linear term (10), we have the uniform estimate 

(16) | (T(g,h),f) | <\g\ L z\h\ N sn\f\ N s„ + |<?U= +2a [\h\v\f\ N .« + Njv^|/U=] . 
For any \a\ < N and N > 3, with the above estimate we have 

\(d a T(g,h),d a f)\ < ||5||^ i ;||h||^ A r-.-r||/ll^^.-r 

+ \\9\\h»L2 +23 [WHh? Ll\\f\\H»N''-i + \\h\\H»N'-v\\f\\H»Ll\ ■ 

This second estimate in Lemma 2 follows easily by integrating (16) over T^, 
using Cauchy-Schwartz and the Sobolev embedding L°°(T^) D H 2 (T^,). The next 
important inequality that we establish is for the linear operator: 

Lemma 3. Consider the linearized Boltzmann operator L = N + K where N is 
defined in (12) and K is defined in (14). We have the uniform inequalities 

(17) |(%3>|<l4„ 

(is) \{Kg,g)\<v\g\l» +c n \g\h, 

where r\ > is any small number and > 0. 

In these estimates there are several things to observe. First of all there are no 
derivatives in the "compact estimate" from (18), which should be contrasted with 
the corresponding estimate in the Landau case [30, Lemma 5] in which the upper 
bound requires the inclusion of derivatives. Further (17) is a simple consequence of 
the main estimate (16). This estimate tells us that the "norm" piece of the linear 
term, given by (Ng, g) , is bounded above by a uniform constant times |g|^r s , T . This 
means that the coercive inequality in the next Lemma 4 is essentially sharp. 

Lemma 4. For the spaces defined in (8) and (13), we have the uniform estimate 

-y + 2s 

The coercive inequality in Lemma 4 and (17) taken together demonstrate that 
the "norm piece" (12) is actually comparable to our designer norm 7V S ' 7 , i.e., 

W. />«!/&..■». 

While the upper bound (17) will have no major use in our arguments, it nevertheless 
is important because it demonstrates that the non-isotropic space -/V s ' 7 naturally 
arises in this near Maxwcllian problem. 

Finally we have a useful, intermediate inequality which is also the key to the 
estimate for K in (18). For the following lemma, we take ej to be a Maxwellian, or 
a Maxwellian times any polynomial p in v of degree at most two: 

Lemma 5. For all positive T),S, we have the following uniform estimate 

[<r(^, cO, />! + |<rCff, /), c,>| < |flrU 1(5 |/I^. 

Functions of the sort as e;, will come up in (78) later on. Of course, Lemma 5 is 
expected to hold whenever e/ is smooth and rapidly decaying, for example, but we 
restrict e; to have the specific form mentioned here because it simplifies the proof 
and the additional generality will be of no use herein. 

In the following we will outline several key new ideas in our proof. 
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2.3. Outline of the article, and overview of our proof. It has been known 
to the experts for some time that the sum total of the inequalities in Section 2.2 
would be sufficient for global existence [32], although crucially the spaces in which 
these inequalities should be proven was unknown. Notice that Lemma 3 implies 

(19) (Lg,g) > |ff|jv«,-r — lower order terms, 

This coercive lower bound inequality is fundamental to global existence. Since the 
the operators T and L are intimately connected, among other consequences, this 
means that if both of (19) and Lemma 2 are simultaneously true, then the Hilbcrt 
space iV s ' 7 satisfying these inequalities is unique. From this point of view the first 
major difficulty which we had to overcome was the identification of the appropriate 
Hilbcrt space (which we have already described). 

Identification of the space N sn . It turned out that the candidate Hilbert space 
iV s ' 7 is a weighted, anisotropic fractional Sobolev space (8) and (12) which corre- 
sponds to fractional differentiation on the paraboloid in R 4 . To estimate this space, 
we find it convenient to use a geometric Littlewood-Paley-type decomposition, in- 
spired by the work of Stein [45]. We do not, however, take a semigroup approach to 
the actual construction of our Littlewood-Paley projections as Stein did. Instead, 
we use the embedding of the paraboloid in R 4 to our advantage. If dfi is the Radon 
measure on R 4 corresponding to surface measure on the paraboloid, our approach 
is to take a rcnormalizcd version of the four-dimensional, Euclidean Littlewood- 
Paley decomposition of the measure gdfi as our non-isotropic, three-dimensional, 
Littlewood-Paley-type decomposition for the function g. Among other benefits, 
this approach automatically allows for a natural extension of the Littlewood-Paley 
projections Pjg and Qjg (from Section 4.1) as smooth functions defined on R 4 in a 
neighborhood of the paraboloid. This allows us to avoid a direct discussion of the 
induced metric on R 3 by phrasing our results in terms of the projections Pjg, Qjg, 
and various Euclidean derivatives of these functions in R 4 instead of R 3 . 

The upper bound inequality. The proof of the main non-linear estimate in 
Lemma 2 is based on a dyadic decomposition of the singularity of the collision 
kernel B in (4) and (5) as well as a Littlewood-Paley-type decomposition of the 
functions h and /. The end result is that one is led to consider a triple sum 

oo oo oo 

E EEl<r*(ff,V>/;>l- 

k— — oo j' — J— 

Here Tfc is the non-linear operator (10) summed over a special dyadic decomposition 
of the singularity, and hj/ , fj are the functions h, f expanded in terms of the 
anisotropic Littlewood-Paley decomposition described just above. Control over the 
sum of the pieces rests on two important observations. First, when considering 
terms for which 2~ fe is large relative to 2~ J and 2~ 3 , a favorable estimate holds 
simply because the support of Bj-Qv — w*|,cos#) is compact and bounded away 
from the singularity at 8 = 0. This is the regime which may be thought of as being 
far from the singularity. Second, when either 2~ J or 2~ J is large relative to 2~ fe , 
i.e., near the singularity, an improvement may be made by exploiting the inherent 
cancellation structure of IV The cost which must be paid in order to use this 
cancellation is that derivatives must fall on either hj' or fj. In this case, with the 
dual formulation (described next) it is always possible to arrange for the derivatives 
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to be placed on the function of our choice. Placing the derivatives on the function 
of largest scale (that is, the function whose index is least) gives some extra decay 
that allows one to sum all the terms by comparison to a geometric series. The 
cancellation structure is not measured in the usual way, we measure cancellations 
using the metric on the paraboloid. 

It should be noted that our analysis allows us to essentially ignore the dependence 
of r(<7, h) on the function g\ this is a great advantage, as it means that one may think 
of the trilinear form (r(<7, h), /) as a family of bilinear forms in h and / parametrized 
by the function g. This observation is essential, since the fully trilinear form falls 
well outside the scope of existing tools in harmonic analysis. 

The dual formulation. A key point of significant technical importance in the 
proof of the upper bound inequality is that we must be able to make estimates for 
(T(g,h),f) which exploit the intrinsic cancellations at the cost of placing deriva- 
tives on any one of the two functions h or / that we choose. If we were not forced 
to consider fractional derivatives, a suitable tool would be integration- by-parts. As 
it stands however, it is necessary to find two different yet analogous representations 
of the trilinear form (T(g, h), /) which clearly relate cancellation to smoothness of 
h and /, respectively. It turns out that placing derivatives on / is fairly straightfor- 
ward to do using existing Carleman-type representations for the bilinear operator 
r. In particular, one may apply a standard pre-post change of variables on the gain 
term Q + to obtain the representation 



which is justified by approximation of B by a sequence of cut-off kernels. Clearly, 
for each fixed g, there is an operator T g such that (T(g, h), f) = (T g f,h), and 
moreover, the formula above can be used to write down an explicit formula for 
T g . To place derivatives on h, on the other hand, it is necessary to find a new 
representation which involves only differences of h' and h, i.e., no differences of g 
or /. To that end, there is a need to compute what we call the "dual formulation," 
since this amounts to writing down a formula for T*. These computations may be 
found in the Appendix; the end result is that 



An interesting consequence of this formula is that the gain term Q + is unchanged 
and only the loss term Q~ differs in these two formulas. These two formulas also 
demonstrate the essentially straightforward dependence on g which we use to apply 
traditionally bilinear methods to the trilinear form. 

The coercive inequality. The key to proving (19), on the other hand, is to show 
the equivalence between (8) and the inner product (Nf,f) from (12). We prove 
equivalent estimates in terms of the Littlewood-Paley projections. This analysis 
consists of two parts. The first is rewriting (13) with a Carleman representation as 



(T(g,h),f) = f dv [ do* [ 



daB g*h{M'J' - M.f), 



JR 3 JR 3 J§ 2 




(20) 
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for an appropriate function K(v, v'), see (60) . If we let d(v, v') denote the Euclidean 
distance in R 4 between the points (v, ^\v\ 2 ) and (?/, ^Iv'l 2 ), a simple pointwise 
estimation of this function K demonstrates that 

K(v,v') > ((«) (v')) 2± ^ ±1 (d(v,v'))- 3 - 2s , 

for a large set of pairs (v, v'), the exact description of which is slightly complicated. 
The second part is to demonstrate that the set of pairs for which this inequality 
holds is large enough to conclude an integral version of this inequality, namely, 

(Nf,f)> [ dv [ dv' ((v)(v')y-^(f'-ff(d(v,v'))- 3 - 2s l d(v , vl)sl . 

This latter argument is accomplished by means of a partition of unity and Fourier 
analysis, the key point being that the expressions 

V J R 3 \v-v'\ 3 + 2s 

are uniformly comparable for all suitable / as Q ranges over the family of nonnega- 
tive, homogeneous functions of degree for which |fi|z,i(§2) > 1 and |f2|£°°(§2) < 1. 

The plan of the rest of the paper is as follows. In Section 3 we will formulate 
the first major physical decomposition of the trilinear form associated with the 
non-linear collision operator (10). With this we prove the main "size and support" 
estimates. We finish this section by formulating the main cancellation inequalities 
using the metric on the paraboloid. 

In Section 4, we develop the anisotropic Littlewood-Paley decomposition which is 
associated to the geometry of the paraboloid. We further prove estimates connecting 
the Littlewood-Paley square functions with our norm (8). 

In Section 5 we prove the key estimate for the trilinear form, Lemma 2. This 
estimate will rely heavily on all of the developments in the previous two sections. 
The "compact estimate" in Lemma 5 will follow shortly from these developments, 
and also the sharp linear upper bounds from Lemma 3. 

The last estimates on the velocity side are contained in Section 6, where it is 
shown crucially that the main norm (8) is comparable to both our anisotropic 
Littlewood-Paley square function and also the space which is generated by the 
linearized operator: (Nf,f) below (12). This involves several ideas, including esti- 
mating a Carleman-type representation and what we call a "Fourier redistribution" 
argument. We further develop useful functional analytic properties of iV s ' 7 . 

In Section 7, we show that all of our singular estimates on the velocity variables 
from the previous sections can be included in the current cut-off theory. Specifically, 
we use the space-time estimates and non-linear energy method that was introduced 
by Guo [30-32] . This works in particular because our new arguments for the velocity 
variables outlined above are morally fully decoupled from the argument to handle 
the space-time aspects of the equation. We further remark that our estimates above 
are, in general, flexible enough to adapt to other modern cut-off methods. 

Lastly, the Appendix contains Carleman-type representations and a derivation 
of a "dual formulation" for the trilinear form (95) that is used in the main text. 
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3. Physical decomposition and related estimates 

In this section we introduce the first major decomposition and prove several 
estimates which will play a central role in establishing the main inequality for the 
non-linear term T from (10) and the norm | • | jv^.t - This first decomposition is a 
decomposition of the singularity of the collision kernel. For various reasons, it turns 
out to be useful to decompose b(cos8) from (4) to regions where 8 w 2~i\v — 
rather than a simpler dyadic decomposition not involving \v — The principal 
benefit for doing so is that this extra factor makes it easier to prove estimates on the 
space L 2 /+2s (M.^ l ) because the weight $>(\v — v*\) from (5) is already present in the 
kernel and the extra weight \v— v* \ 2s falls out automatically from our decomposition. 

The estimates to be proved fall into two main categories: the first are various 
L 2 - and weighted L 2 -inequalities which follow directly from the size and support 
conditions on our decomposed pieces (such estimates are typically called "trivial" 
estimates). The second type of estimate will assume some sort of smoothness and 
obtain better estimates than the "trivial" estimates by exploiting the cancellation 
structure of the non-linear term T from (10). It is already worth stating at this point 
that the particular smoothness assumptions we make are dictated by the problem 
and will specifically be somewhat unusual; in particular, they will not correspond 
to the usual, Euclidean Sobolev spaces on R 3 . 

3.1. Dyadic decomposition of the singularity. Let {Xj}j^L-oo be a partition 
of unity on (0,oo) such that |xj |l°° < 1 and Xj ' 1S supported on [2 _J_1 , 2~ J ']. For 
each j , let 

Note that 

i /i2 \ v ~ v *\ 2 fi / v - v* \\ 2 . 2 6 



, a ) = \v — v*\ sin 



Hence, the condition \v — v'\ « 2~- 7 is equivalent to the condition that the angle 
between a and Z-v*\ 1S com P ara ble to 2 _ - 7 |v — With this partition, we define 



TU9,h,f) 
Ti(g,hJ) 



/ dv dv* I da Bj(v — v*, a) gJiM'^f 
Jr 3 Jm. 3 Js 2 

dv dv* / da Bj(v — v*,a) g*hM t f. 

JR 3 JR3 J§2 



It turns out that we will also need to express the collision operator (10) using its 
"dual formulation." With the variant of Carlcman's representation coming from 
Proposition 10 and the notation M* = M{v + v„ — v') ( = M A A /,' on E% t ), we have 
the following alternative representation for Ti_ as well as the definition of a third 
trilinear operator Ti (based on the calculation (95) with, recall, v^ = v + v* — v'): 

T{(g, h,f) = if dv' [ dv t [ dn v Bj(v - v.,2v' - v -v,) , 
Jr3 Jr3 Je*[ F - v *\ \ v ~ v *\ 

Ti(g,hJ)^U[ dv' [ dvj dn v B 3 |^ = v *\ ^— ^ g*f'M t h', 
Jr 3 Jr 3 Je?,' $(u-u*J \v-v*\ 4 
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where above we make the slight abuse of notation that 

v — i>* 1v' — v — v* 



Bj = Bj(v-v*,2v'-v-v*) = <&(\v-v*\)b 



\v — ' \2v' — v — v*\ 



In these integrals above djr v is Lebesgue measure on the two-dimensional plane E% 
passing through v' with normal v' — , and v is the variable of integration. When 
f,g, h <G 6>(R 3 ), the pre-post collisional change of variables, the dual representation 
(95) from the appendix, and the previous calculations guarantee that 

oo 

(T(g,h)J) = £ {Ti(g,h,f)-TL(g,h,f)} 

j=-oo 

OO 

= {Tl(g,h,f)-Ti(g,h,f)}. 

j=— oo 

These will be the general quantities that we estimate in the following sections. The 
first step is to estimate each of Ti. , T°_ , and Tl using only the known constraints 
on the size and support of Bj. 

3.2. "Trivial" analysis of the decomposed pieces. Under the condition that 
7 + 2s > 0, the following basic inequality holds uniformly for all 77, 5 > 0: 



(21) 



We will use this inequality in all of our estimates and refer to it as the (77, <5)- 
incquality. We begin with the following: 

Proposition 1. For any integer j and any rj, 8 > 0, we have the uniform estimates 



(22) 
(23) 



T 3 l(g,h,f) <2 2sj \g\5 ir> \h\ v , s \f\ L * +2s 
TL(g,h,f) <2 2s ^g\ s , v \h\ Ll \f\ v , s . 



Moreover, uniformly for any any integer j we have 



(24) 



TL{g,h,f) <2 2 ^\g\ L ,\h\ Ll \f\ 1 



Proof. Given the size estimates for b(cos9) in (4) and the support of Xj> clearly 

de e- 1 - 23 < 2 2s i (v - i7*) 7+2s . 

(Note that this inequality holds true when $(w) equals either |w| 7 or (i>) 7 because 
of the assumption 7 + 2s > 0). Thus 



(25) 



s 2 



da Bj < <l>(\v-v*\) 



Ti(g,h,f) <2 2s > / dv dv*M*(v-v*) y+2s \g*\\hf\. 



Taking a geometric mean of the (77, i5)-inequality for (v — v*) y+ s , i.e., (21), and the 
inequality ill* (v — v*)~ f+2s < (v)~ f+2s gives that 

mi (v - v*y +2s < (( v y +2s ( v ( v y +2s + s- 1 ) (s (v*y +2s + r,- 1 



With Cauchy- Schwartz on the 7j* and v integrals, for any j and any 77 > 0, 

Ti(g,h,f) <2 2 ^\ g \ s , 1 \h\^\f\^ +2s - 
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Since T J _(g,h, f) is symmetric in h and /, (22) must imply (23). As for (24), note 
that the inequality M^ 2 (v — v*) 1+2s < (w) 7+2ji implies (24) by exactly the same 
reasoning used in the first two inequalities. □ 



Proposition 2. The inequalities below are uniform for any integer j and r/,5 > 0: 

(26) \Tl(g,h,f)\ <2 2 v\g\ s , v \h\ v , s \f\ L 2 4 

(27) \Ti(g,hJ)\<2 2 ^\g\ Str ,\h\ L 2 +2 
Moreover, the following inequality also holds uniformly for any integer j : 



\fks- 



(28) 



\Ti(g,h,f)\<2 2 ^\g\ L2 \h\ 1 



.l/l; 



Proof. As in the previous Proposition, the key to both inequalities is the symmetry 
between h and / coupled with two applications of Cauchy-Schwartz. The difference 
is that, this time, the Carleman representation will be used and the main integrals 
of will be over v* and v' . In this case, the quantity of interest is 



c?7r„ b 



\v'-v*\ 2 - \v-v'\ 2 \ \v'-v*\ 



\v' — u*| 2 + \v — v'\ 2 J \v — u*| 



- Xj (\v-v\). 



The support condition yields \v — v'\ w 2 3 . Moreover, since b(cos8) vanishes for 
9 e [7r/2,7r], we have \v' — u*| > \v' — v\. Consequently, the condition (4) gives 



\v' — v*\ — \v — v'\ 



/|2 



< 



\v - v'\' 



^ ^ \ \v' — u*| 2 + \v — v'\ 2 

Thus, the integral is bounded by a uniform constant times 
\v'-v*\ 2+2s 



dir„ 



E y' 



1 ,/|2+2s 



\v'-v4- 2 Xj (\v-v'\)<2 2s i\v'-v*\ 



2 s 



As a result 



\T>(g,h,f)\ 



< 



dv' / dv*M*\g*ti f'\2 2sj (v' - v*) 



7+2s 



This leads directly to (26) and (27) in the same way that (22) and (23) were 
obtained. Similarly, (28) follows from M} /2 (v* - v') 1+2s < (z/) 7+2 '\ □ 

Proposition 3. For any integer j and 5,r) > 0, we uniformly have the estimate 



(30) 



Tl(g,h,f) 



9\s, n \ h \vAf\L^ 



Moreover, for any j > we have the uniform estimate 



(31) 



Tl(g,h,f) <2 2s i\g\ L2 \h\ 1 



.\f\i 



Proof. The proof of (30) follows by interpolation. First, we will assume that / G 
L°°. Noting that M' t < 1, the integral of Bj over § 2 may be estimated with (25) 
to give that 



Tl(g,h,f) < 2 2 »\f\ L ~ f dv f dv.\ 

JS. 3 JR3 



(v ~ v*) 



7+2s 
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Using the (77, <5)-inequality for (v — v. 
Ti(g,hJ)\<2 2s i\f\ L < 

(32) 



For both g and h, this is exactly the .L -type estimate which will give the | ■ \s tV or 
the I • \ri t s norm, respectively, after interpolation with L°° . 

For the next estimate we use g,h G L°°. In this case, the pre-post collisional 
change of variables yields 



(inequality (21)), it follows that 
dv.\g.\ (5 (v*)~< +2s + rj- 1 ) 

dv\h\{r) (w) 7+2s +5- 1 ] 



Ti(g,h,f) 



dv / dv* / doBjM^g'Jh! . 



Then estimating the integral of Bj as in (25) grants 



Ti(g,h,f) < 2 2s i\g\ L ~\h\ L ~ / dv dv*M*(v-v*) 



\7+2s 



l/l- 



After exploiting the inequality M* (v — 

y+2s < M l/2 ^7+2^ the integral 

over 

may be bounded above to yield 



(33) 



rjteA/) < 2 2sj \g\L~\h\ L ~ f dv \f\(v) 



7+2s 



Interpolating (32) and (33) gives (30) (the Riesz-Thorin interpolation theorem with 
weights found in Stein and Weiss [46] suffices here; since g and h are in the same 
space, we do not even need a multilinear interpolation theorem since we can instead 
treat g*h as a function on K 3 x R 3 which is either in a weighted L 1 space or in L°°). 

Regarding (31), the proof is essentially unchanged. The only difference is, when 
/ is taken in i°° and g, h are taken in L l , one may improve the inequality < 1 
to Ml < M 4 1_e for any fixed e > because \v* — v'J = \v — v'\ < 1 when j > 0. 
Setting 77 = 6 = 1 gives the result, since any wcig ht (v*) 13 Ml~ e is bounded. □ 

Proposition 4. The following inequality holds uniformly in j and 5,rj > 0: 



(34) 



TUg,h,f) 



<2^\g\ s , v \h\ L . \f\ n , 5 . 



Proof. The proof proceeds along the same lines as above, the difference is that this 
time we use the Carleman representation for T|. In this case, as in (29), we have 

2v' — v — u* 



Bj ?J — V* , 



2v' — v — vA 



\v — v l 2+2s 

<*(«-«^l«-«'l) 



\v — v'\ 



/|2+2s 



< - v*)Xj{\v - V 



1 1 \ \v — vA, 



2+2s 



\v — v'\ 



/\2+2s 



The inequalities above hold because \v — v* | 2 w \v' — u* | 2 and as in the analysis of 
(29) we know that \v' - vA\ 2 > \v - v'\ 2 . Thus 



dn v-, n— 7 

E v> \v-v*\\v' 



< - V*)\v' - 7J*| 



2 s 



6?7T, 



ex' 



\v — V 



/|2+2s 



< 2 2s] (v' 



\l+2s 
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As in the previous cases, this leads to the inequality 



Tl(g,h,f) <2 2 »\h\L-([ 



dv*\g*\(8 (v*) 7+2 ' s + 77 1 

dv'ifKvivy^ + s- 1 ) 



In a similar manner, after a pre-post change of variables, we obtain 

Ti(g,h,f) <2^\g\ L ~\f\ L ~ [ dv' \h'\(v') 1+2s ■ 

We finish the proof of (34) by interpolation in exactly the same manner that the 
corresponding inequality (30) was obtained in the previous proposition. □ 



3.3. Cancellations. In this section, we seek to establish estimates for the differ- 
ences Tj? — T_ and Ti — T*. We wish the estimates to have good dependence on 
k (in particular, we would like the norm to be a negative power of 2 k ), but this 
improved norm will be paid for by assuming differentiability of one of the functions 
h or /. The key obstacle to overcome in making these estimates is that the mag- 
nitude of the gradients of h and / must be measured in some nonisotropic way; 
this is a point of fundamental importance, as the scaling is imposed upon us by the 
structure of the "norm piece" (Nf, f) . 

The scaling dictated by the problem is that of the paraboloid: namely, that the 
function f(v) should be thought of as the restriction of some function F of four 
variables to the paraboloid (v, ^\v\ 2 ). Consequently, the correct metric to use in 
measuring the length of vectors in R 3 will be the metric on the paraboloid in R 4 
induced by the four-dimensional Euclidean metric. To simplify the calculations, we 
will work directly with the function F rather than / and take its four-dimensional 
derivatives in the usual Euclidean metric. This will be sufficient for our purposes 
since our Littlewood-Paley-type decomposition will give us a natural way to extend 
the projections Qjf into four dimensions while preserving the relevant differentia- 
bility properties of the three-dimensional restriction to the paraboloid. 

To begin, it is necessary to find a suitable formula relating differences of F at 
nearby points on the paraboloid to the various derivatives of F as a function of 
four variables. To this end, fix any two v, v' <E R 3 , and consider 7 : [0, 1] — > R 3 and 
7 : [0, 1] -> R 4 given by 

7 (6») = f 0v' + (1 - 0)v, and j[6) = (ov' + (1 - 6)v, i \6v' + (1 - 9)vf 

Note that 7 lies in the paraboloid {(«i, . . . , U4) € R 4 U4 = \{ v \ + ■ ■ ■ v l) } , an d 
that 7(0) = v and 7(1) = v' . Elementary calculations show that 



d^y d^^i 

-J: = (y'-v,{>y(p),v , -v)), and = (0, \v' - V 



2\ 



Now we use the standard trick of writing the difference of F at two different points 
in terms of an integral of a derivative (in this case the integral is along the path 7) : 



/" • (v^XtW) 
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where the dot product on the right-hand side is the usual Euclidean inner-product 
on M 4 and V4 is the four-dimensional gradient of F. For convenience we define 



\V 4 \ k F( Vl ,..., Vi)= sup (f- V 4 rF(ui,...,u 4 ) 



l«l<i 



0,1,2, 



where (el 4 and |£| is the usual Euclidean length. In particular, note that we 
have defined |V 4 |°-F = \F\. 

If v and v' are related by the collision geometry (2), then (v — v' , v' — v m ) = 0, 
which yields that 

(j{6),v' -v) = (v*,v'-v) - (1 -6)\v-v'\ 2 . 
Thus, whenever \v — v'\ < 1, which is the case of interest, we have 



dO 



< \v-i/\(v m ) 



Indeed, throughout this section we suppose that \v 



< 1 since this is the 



situation where our cancellation inequalities will be used. In particular, we have 
the following inequality for differences related by the collisional geometry: 



(36) 



F[v', 



,'|2 



F[v, 



<(v*)\v-v'\ / d6 |V 4 |F( 7 (#)). 



Furthermore by subtracting the linear term from both sides of (35) and using the 
integration trick iteratively on the integrand of the integral already appearing on 
the right-hand side we obtain that 



(37) 



F v' 



„'|2 



- J(0)-V 4 i^) 
1 



< 



(v.) 2 \v-vf / dO [|V 4 |F(7(fl)) + |V 4 | 2 F(7(0))] 



We note that, by symmetry, the same result holds when the roles of v and v' arc 
reversed (which only changes the curve 7 by reversing the parametrization: j(9) 
becomes 7(1 — 6)). We will use these two basic cancellation inequalities to prove 
the following cancellation estimates for the trilinear form: 

Proposition 5. Suppose f is a Schwartz function on R 3 given by the restriction 
of some Schwartz function F on R 4 to the paraboloid (v, ^\v\ 2 ). For each k, let 
|V 4 |' C / be the restriction of |V 4 | fc i 7 ' to the same paraboloid. Then, for any j > 0, 



(38) 



pi-TL)(g,h,f) <2( 2s - 1 ^| 5 | L2 |/ i | i 



Ei V4 i"/ 



fc=0 



If s > \, then it also holds uniformly for all j > that 



(39) 



(Tl-Ti)(g,hJ) <2( 2s -^| 5 | L2 |/ l | i 



Eiv 4 iv 



fe=0 



Proof. Write the difference M'J' ~ M m f = M*M ((M -1 /)' - (M -1 /)). If / ex- 
tends to M 4 , then (M _1 /)(^i, • • • ,^4) = e Vi/2 f(vi, v 4 ). With this extension 

(40) |V 4 |(M-7) < M^IVil / + Af- 1 |V 4 |/. 
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To prove (38), we first observe from (36) that for any fixed e > we have 



\{Tl-Tl){g,h,f)\<2-n dd dv dv* d*B j Ml- c \g.\\h\f(0v'+(l-9)v), 

JO JR 3 JR 3 JS 2 

where / = |V 4 | / + |V 4 |/. The loss of e arises from the factor of (v*) in (36), 
which also accounts for the factor of 2~- 7 . Additionally we have used the inequality 
M'^M' 1 = MM'' 1 < e 2 ~ 3 \ v >\ which follows from 

(41) \v'\ 2 - \v\ 2 = -2 (v', v-v')-\v- v'\ 2 = -2 (v„v - v') - \v - v'\ 2 . 
By Cauchy-Schwartz, it suffices to prove the estimates 

i 

(42) ( f d6 f dv [ dv, f dtTB j Mi- e \g.\ 2 \h\ 2 Y < r*\g\v\h\ L * , 

\J0 JR 3 JR 3 JS 2 ) 7 + 

and 

(43) ( [ d6 [ dv f dv, f daB 3 Ml^\~f(9v' + {1 - 6)v)\ 2 Y <2^\f\ L , . 

\J0 JR 3 JR 3 J$ 2 / ~' + 

The former inequality, (42), follows from (25) as before. The latter uniform bound, 
(43), follows after the well-known change of variables u = dv' + (1 — 6)v, which 
changes v to u. With (2), we see (with Sij the usual Kronecker delta) that 

dVj v ' 3 dVj V 2 / 2 
with the unit vector k = (v — — Thus the Jacobian is 

dm ( e\ 2 r / e\ e „ . 



dvj 

Since b((k,a)) = when (k,a) < from (4), and 9 G [0,1], it follows that the 
Jacobian of this change is bounded from below on the support of the integral (43). 
But after this change of variable the old pole k = (v — v*)/\v — v*| moves with the 
angle a. However it is easy to check that, when one takes k = (u — v*)/\u — i>*|, 
1 — (fc, a) ss 1 — (k,a) , meaning that the angle to the pole is comparable to the 
angle to k (which does not vary with a). Thus the estimate analogous to (25) will 
continue to hold after the change of variables, giving precisely the estimate in (43). 

The proof of the inequality (39) proceeds in exactly the same fashion, using (37) 
instead of (36). In this case, similar to (40), pointwise everywhere in R 3 we have 
|V 4 | 2 (M" 1 /) < il/- 1 |V4| / + M- 1 |V 4 |/ + M- 1 |V4| 2 /. By subtracting off 

§(0)-V 4 FH 

we easily obtain (39), at the price of still needing to estimate this term alone. Here 
the extension of A/ -1 / is once again inserted in the place of F in (37). 

Notice that (0) is linear in v' — v and has no other dependence on v'; af- 
ter multiplying it by Bj and integrating with respect to er, the symmetry of Bj 
with respect to a around the direction |"Z"*| forces all components of this integral 
to vanish except the component in the symmetry direction; in other words, one 
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may replace v — v' with ^„"* | (v — v' , \ v v Jl* | \ in the expression for ^f(O). Since 

— 'i 2 

: 0, the vector further reduces to njE^fi \ v - v \ ■ Simply observing 



v — v* \v — v'\ 



\v — u* \v — v*\ 



< 2- 2j \v-v*\~ 1 , 



allows one to employ the same methods as above (easier in this case) to estimate 



this term for the projection of ^(0) onto the first three of the four coordinate 
directions. To be precise, one must bound the following integral 

(44) 2~ 2j [ dv f dv* f daBjMl-^Wh^V^f +\V i \ 1 f)\v-v*\-\ 

JR3 JR3 JS 2 

Here, once again, we absorb any powers of by M~ e . The estimation of this 
integral proceeds exactly as was done for (24), the only difference being the extra 
factor \v — In particular, (44) is bounded above by 

2( 2s ~ 2 ^ j dv f ^$( W -^)| U -^| 2s - 1 Mi- e |^||ft|(|V4| / + |V 4 | 1 /)- 

When \v — w*| > 1, \v — w»| 2s_1 < (v — v*) 2 ' s , and we obtain the upper bound 

k=0 

which has a weight on h and / which is even better than desired. As for the 
remaining piece where \v — v„\ < 1, we have the upper bound 



dv r Mt- e \g*Mv - v*)\v - v.r- 1 !^,^ < \Mf e g*\ L 2, 

uniformly in v by Cauchy- Schwartz as long as 7 + 2s — 1 > ^ so that the L 2 - 
norm of the weight <£(i> — u*)|v — v *\ 2s ~ 1 l\v-v t \<i in the variable v t is uniformly 
bounded above as a function of v. Integrating in v and applying Cauchy-Schwartz 
to separate h and / gives the desired inequality on this piece as well. 

The fourth coordinate direction of S(0) is given by (v,v' — v), which reduces to 



v — v* \v — v'\' 
\v — vJ \v — vJ 



< 



\v - v'\' 



Therefore this term can also be handled in the same way to the previous cases. □ 

Proposition 6. As in the previous proposition, suppose h is a Schwartz function 
on R 3 which is given by the restriction of some Schwartz function in R to the 
paraboloid (v, ^\v\ 2 ) and define IV4 analogously. For any j > 0, the inequality 



(45) \TUg,hJ)-Tl{g,h,f)\<2^-^\g\ L2 \f\ 1 



fc=0 



Moreover, if s > h, then uniformly for j > we have 
(46) \TUg,h,f) -Ti(g,h,f)\ <2 2 ^\g\ L2 \f\ Ll _ 



J22- kl - ekj \V±\ k h 



k=0 



where > and 2s — k — e& < for k = 0, 1, 2. 



GLOBAL SOLUTIONS OF THE BOLTZMANN EQUATION WITHOUT CUTOFF 



2.3 



Proof. This proof follows in the same pattern that is by now well-established. The 
new feature in this case is that the pointwisc difference to examine is 

A = 4 {Mh${v - v*)\v - v»| 3 - M'ti$(v' - «*)|«' - w*| 3 } ■ 

As before, we study the function of 9 given by 

(M/i)(7(0))$( 7 (0)-^)| 7 (0)-^| 3 , 

since A is a constant times the difference in the values of this function at 9 = 
and 9 = 1 (note that we also employ the extension of Mh to M 4 and make the de 
facto extension of the factor $(i> — v*)\v — w*| 3 to M 4 by assuming the new function 
is constant in the fourth variable). In terms of A, our operator may be written as 



We use (36) and (37) again to estimate A. The additional fact required here is that 

|v 4 | fe ($(« - «.)!« - «*l 3 ) < |« - «*!"*(*(« - - u *l 3 )> 

which simply comes from differentiating with respect to v (since the extension is 
taken to be constant in the fourth direction, the gradient V4 reduces to the usual 
three-dimensional gradient). Applying (36) gives the estimate 

d9\{Mh){ 1 {9))\ 



|A| 



<E>(v — v*)\v — w*| 3 



1 

d9 (M|V 4 | /i + M|V4| 1 /i)(7(^)). 







Again, \v — v'\ < \v — v* \ and \v — v'\ < 2 J , so as in the previous Proposition, we 
may estimate the difference T J + — T* by a sum of three terms: 



I c = I d9 I dv' I dv* I dn v -. rBj(v — u*. 2v' — v — «*) 



\v — v* \\v — V* 



7= f 2" J ' / d0 



df' / dw* / dir v Bj(v — v*,2v' — v — v*) 
Jr 3 Je^' 

Ml- e \g*f 



v*\\v 



Va\ h(9v' + (l-6)v) 



I = 2- j 1 1 d9 



dv' / dv* / dn v Bj(v — v*, 2v' — v — v*) 
Jr 3 Je?.' 



x 



Mt e \9*f\\V4\h( 7 (6)) 



v 



In each case, the extra factors of (v*) and M(M') -1 are absorbed into a single 
factor M~ e , possible because of (41) and the fact that j > 0. Each of the terms 
II and III is completely analogous to a corresponding quantity which arose in the 
previous proposition. Splitting M\~ e g i( j'h = (All~ e ^ 2 g*f')(M^ 2 h) allows one to 
employ Cauchy-Schwartz just as was done for (42) and (43) to separate the estimate 
into one integral involving only g and / (which may be estimated using exactly the 
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same argument that handles (28)) and one integral involving only h. The integral 
involving only h is handled by changing from the Carleman representation back to 
the sigma representation. Since $(?/ — v*) ~ $>(v — v*), the desired result follows 
directly from the same argument used to estimate (43). As for term J, the two 
relevant inequalities to establish are 

\I\<2^-^\g\ L2 \f( V y +2s \ L1 \h\ L ^ 

\i\<2 (aa - 1)j \g\v\f\L-\h(vy +2 '\ L i. 

The first is completely analogous to the extra estimate in the previous proposition 
to handle s > \ and follows as long as 7 + 2s — 1 > just as before. The second 
may be transformed into exactly the same form after reverting from the Carleman 
representation to the sigma representation and using the same change of variables 
from the previous proposition as well. Finally, these two results are interpolated to 
give the desired L 2 inequality. 

Regarding the proof of (46), once again, we split the estimate into two parts: 
one involving the linear correction term and one involving the integral of the second 
derivative. This latter term may be estimated by a sum of expressions of the form 



r r r U,_^|2 



d6 / dv' I dv* / dir v rrBj(v — v*, 2v' — v — v*) 

Jk3 Jw* Jev! \v-v,\* 



Ml- e \g*f'\ 
' — vJ\v' — vJ 



\V 4 \ e h{6v' + (l-6)v), 



where k + £ < 2. Now 



\v — v*\" \v — v*\* \v — v*^ 

for any p <G [0, k}. Choose p so that 7 + 2s — p > —I and 2s — 2 + k — p < 0; 
note that this is always possible when 7 > — |. Using this estimate and breaking 
into regions where \v — v*\ < 1 and \v — u*| > 1 allows one to estimate the piece of 
T J + — Ti governed by the right-hand side of (37) by a sum 
2 

E 22s3 >U 2 i/i^ +2s 2-«- e "iiv 4 i^u3 +2s , 

where eg is some positive number (namely, 2 — £ — k + p) such that 2s — i — e; < 0. 
We remark that improvements can be made to this estimate by instead taking the 
norm of h when I = and k = 2 and, after a Sobolev embedding, estimating 
isotropic derivatives of our anisotropic Littlcwood-Palcy decomposition. With that, 
this estimate can be made to work for all 7 + 2s > —1 as well. 

Last, but not least, is the analysis of the linear correction on the left-hand side 
of (37). Without loss of generality, we may instead choose to use the linear term 

^(1)-V 4 F(7(1)), 

in (37) instead of evaluating at 9 = 0. Again, we exploit the symmetry of the kernel 
Bj in the plane E% around the point v'. To simplify matters, the equality 

|=(1) • V 4 F(7(1)) = M'*(t/ - - t'*| 3 ^(l) ' (W - \e 4 h>\ , 



GLOBAL SOLUTIONS OF THE BOLTZMANN EQUATION WITHOUT CUTOFF 



25 



does not involve derivatives of $(|t/ — — v*| 3 by virtue of the fact that 

^\ 1 (6)~v*\ 2 = 2(v'-v, 1 (8)-v*) = 0, 

when 9 = 1 because (v — v', v' — v*) =0 (here is the unit vector pointing in the 
fourth direction in R 4 ). We are therefore left to estimate the integral 



dv' / dv* I d-K v Bj (v — , 2v' — v — v* ) - 
Jts? Jev' 



This integrand still has the property that as v varies on circles of constant distance 
to v , the entire integrand is constant except for (1). If we write ^jf (1) as a sum 
of two vectors, one lying in the span of the first three directions and the second 
pointing in the fourth direction, it follows that we may replace the former vector 
by its projection onto the direction determined by v' — v„. But since the original 
vector points in the direction v — v', the projection vanishes. In other words, only 
the projection of -js (1) pointing in the fourth direction remains. Its magnitude in 
this direction is exactly (v 1 , v' — v), the corresponding integral of which over v also 
vanishes by symmetry. Thus, the linear term integrates to zero in this case. □ 

Finally, let us observe that for j > 0, we have the following uniform inequalities 

2 



(47) \Ti{g,hJ)-TL{g,h,f)\<2^-^\g\ 5 , n \h\ v , s 



fc=0 



(48) \TUg,h,f) -Ti (g, h, f)\ < 2 2 ^\g\ 5 , n \f\ n , s 



2 



J22~ k] - ekj \\/ 4 \ k h 



k=0 



These estimates follow immediately from the work above after taking into account 
(15). This observation will be necessary to establish favorable estimates for the 
"compact piece" given by (T(g, M), g) as well as for the term (T(g, g), M) . 

4. The anisotropic Littlewood-Paley decomposition 

In this section we introduce the Littlewood-Paley decomposition. Rather than 
the standard Littlewood-Paley decomposition on K 3 , we will instead use a decom- 
position which is implicitly adapted to the induced Laplacian on the paraboloid 
(v, ^H 2 ) C R 4 . The main reason for doing so is that an analysis of the norm piece 
(Ng,g) from (12) shows an inherent anisotropy in the directions of differentiation. 
(The analogous problem for the Landau equation is anisotropic [30] , and the sharp 
norm for the non cut-off problem was conjectured to be anisotropic in [42].) Rather 
than work directly on the paraboloid, though, it turns out to be somewhat simpler 
to think of the Littlewood-Paley decomposition we use as being a 4-dimensional 
Euclidean decomposition restricted to the paraboloid. In the process, we will effec- 
tively construct a corresponding extension operator for bandlimitcd functions (i.e., 
functions with dyadically localized frequency support) on the paraboloid into ban- 
dlimitcd functions on R 4 in the neighborhood of the paraboloid (and the extension 
operator satisfies favorable Sobolev-type inequalities). This trick will allow us to 
use the estimates from the previous section without any serious concern for the 
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deeper geometric aspects of our nonisotropic construction (which contrasts with 
the approach of Klainerman and Rodnianski [36]). 

4.1. Associated calculations and definitions. Throughout the remainder of 
this section, we will use the variables v and v' to refer to independent points in 
R 3 , meaning that we will not assume in this section that they are related by the 
collision geometry. The reason we choose to use these variable names is that they 
give a hint about where the Littlewood-Paley projections will be later applied in 
situations which involve the collision geometry explicitly. 

First we develop some tools for calculus on this paraboloid. For v € R 3 , let 
v = f (v, i | z; | 2 ) <G R 4 . Perhaps the most useful such tool is the following: for any 
»eR 3 , let t v : R 3 R 3 and t v : E 3 R 4 be given by 

t v (u) =' u — (1 — (u)" 1 ) (v, u) \v\~ 2 v, and t v u = f (t v u, (v)" 1 (v, u)). 

These mappings should be thought of as sending the hyperplane U4 = to the 
hyperplane tangent to the paraboloid (v, ^H 2 ) at the point v. It's routine to check 
that (v,t v u) = (v) 1 (v,u) and |t„«| 2 = |u| 2 — (v) 2 (v,u) 2 , which implies 

Ir^ul 2 = \t v u\ 2 + (v, t v u) 2 = |ii| 2 , 

meaning that r v is an isometry from one hyperplane to the other. Moreover, it is 
easy to check that 

!| ,2 

V + T V U = V + t v u + - \u\ e 4 , 

where e 4 = f (0, 0, 0, 1). The last term will be thought of as a perturbation, and is 
the basis for all the analysis that follows: 

Proposition 7. Suppose ip is any fixed, smooth junction supported on the unit ball 
in K 4 . For any j > 0, the expansion 

(49) Lp (2 j ( v + 2~ j T v u - v)) = tp(z v u) + 2-*- 1 tp i (r v u) + E s (u, v), 

holds, where ipi = and \Ej(u,v)\ < 2~ 2 i and is supported on a set \u\ < 2. 
Moreover, if \u\ < 1, then 



P-2 



(v + t v u)' - (vf - 2-i p (v)- 1 (v, u) (vf~ 2 
uniformly in v, u, and j > for any fixed /3 £ TSL. 

Proof. Note that (w)' 3 = (1 + 2(u) 4 )' 3 / 2 . With this observation and the equality 

|w| 2 e 4 , 



v + 2- j r v u = v + 2~ j T„u + 2- 2j - 1 \- A2 - 



both inequalities immediately follow from Taylor's theorem with remainder. □ 
These expansions may be utilized to estimate a variety of integrals of the form 

dv' 2 3j (p(2 j (i/_-v)) (v'f = (vy 1 j du ip{2 3 { v + 2- j t v u -v)) (v + 2- j T v u) P . 

Jr 3 

The right-hand side follows by the change of variables v' n- v + 2-3t v u as suggested 
by the proposition above. If it is assumed that (p is radial, then ip^ is odd, meaning 
that lP4:(t v u) is an odd function of u (and, by assumption, <p{r v u) is an even function 
of u). Likewise (v)^ is (trivially) an even function of u, while (v,u) (w)^~ 3 is odd. 
Thus if both factors in the right-hand side of the above equality are expanded as 
a function of u by means of the previous proposition, there are a number of terms 
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which automatically cancel (namely, the product of the zeroth-order term for one 
and the first-order term for the other). Thus the result is that 



(50) 



dv'2 3j <p(2 j {v-y!_)) (v'f - (v) 



duip(z v u) 



< 2~ 2j 



v/3-1 



Moreover, if tp is assumed to be radial, it is also true that J dutp(r v u) is a constant 
(since t v is an isometry of M. 3 with some hyperplane in R 4 passing through the 
origin). It will be assumed that cp is chosen so that this constant equals 1. 
We now define the following anisotropic Littlewood-Paley projections: 



Pjf(v) 



f(v')2 3j ip(2 j (v-v!)) («') dv', j > 0, 



Given the calculations above, if / is a Schwartz function, then clearly Pjf(v) 
as j — > oo by the normalization condition on tp and the estimate (50). That same 
estimate proves the fundamental fact that 



(51) 



dv{P 3 f(v)Y (v 



< 



dv(f(v)) 2 (v) 



uniformly in j > for any fixed j3 € K. and any p £ [l,oo) (as a consequence of 
Schur's test for integral operators). 

4.2. Square function estimates related to the norm (8). The next step is 
to establish a favorable inequality relating the associated square functions, i.e., 
{^2j2 2s: >(Qjf) 2 {v)} 1 l 2 , to an integral involving a squared difference of the form 
found in (8). To this end, notice that 

\fdv(dv'( dz(f(v) f(v')) 2 qj(z,v)qj{z,v!_) (v) (v>) (zf 



K 3 



R 3 



dv(Q 3 f(v)) 2 (v) 



dvQ^n^Q^iv) (v) 



with qj(z,v) = [2 sj (p(2^(v-z)) - 2 3: >- 3 (p(2 : >- 1 (v - z))] . To see this, expand the 
square (f(v) — f(v')) 2 and exploit the symmetry of the integral in v and v'. Next, 
recall that |Q J (l)( , i;)| < 2~ 2 i by virtue of (50) and the triangle inequality. But 



dv\Q j (l)(v)\\Q j (f)(v)\(vf<2- 2 i I dv'(f(v')) 2 {v'^ 



R 3 



where the last inequality follows from the pointwise estimate for Qj(l) and the 
boundedness of Qj on L 1 (i.e., the inequality (51) for p = 1). Consequently, for 
any j > 0, we have the following estimate: 



dv / dv' / dz(f(v) - f(v')) 2 qj {z,v)qM,v!_) (v) (v') (z)' 



dv{Q 3 f{v)) 2 (v) 



R 3 



< 2- 2 > / dv\f(v)\ 2 (vf 



This inequality allows us to compare the weighted norm of Qjf to an integral 
involving the squared difference (f(v) — f(v')) 2 , the right-hand side of the inequality 
above being the lower order error in this comparison. 
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The next step is to estimate the integral of (z)^ qj(v — z)qj(v — z) with respect 
to z. By the same reasoning that lead to (50) (namely, by making the change of 
variables zh>« + 2~ j t v w), it must be the case that the inequality 

dz(zf\ qj (z-v)\<(vf'\ 

holds uniformly for j > (in fact, this follows directly from (50) and the triangle 
inequality). Moreover, clearly \qj(z_— v)\ iS 2 3j as well. Thus, for any j > 0, 



dz (zf qj (v - z)q 3 {z-y^j < 2 3j (vf' 1 , 
and the integral is supported on some set \v — %/_ | < 2" . Consequently 
f> 2 ^ / dv(Q 3 f{v)f{v) fi 



3=0 



4 



where xiv. — vL) is some bounded, nonnegative function supported on a ball in 
given by \v — ?/| < 1. Note that, if the scaling function (p is suitably rescaled on R 4 
(i.e., (p(v) is replaced with (p(cv) where c is the implicit constant in \v— j/| < 1), then 
X may be taken to be the characteristic function of the unit ball in K 4 . Furthermore, 
\v— 1/| < 1 implies (v) ~ (v 1 ) so the powers of v and z/ may be redistributed at will. 
In particular, then, it must hold that the exponentially weighted, squared sum of 
Littlcwood-Paley Qj's is dominated by the square of our fundamental norm (8): 



£2 2 ^ / dv(Q 3 f(v)f(vy +2s <\f\ 



Likewise, since each Qj has a natural extension to M 4 (obtained by replacing v by 
an arbitrary 4- vector in the definitions of Pj and Qj ) , it is natural to ask a similar 
question about the 4-derivatives | V 4 1 fc Q ^ . In this case, it is easy to see that the 
same estimates must hold (as the derivatives simply fall on 93), except that a factor 
of 2 J is introduced for every derivative. In particular, then, we must also have 

(52) ^ 2{s - k) ° / dv(\V 4 \ k Qjf) 2 (v) (vy +2s < l/lk,, k = 0,1,2, 

uniformly in j > as well. These inequalities will be fundamentally important in 
the next section, where all the estimates so far are combined to yield the stated 
upper bound of the trilinear form (T(g,h),f). 



5. Upper bounds for the trilinear form 



In this section, we establish Lemma 2 for the nonlinear term T as well as Lemma 
3 which provide estimates related to the splitting L = N + K of the linear term. 
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5.1. The main upper bound inequality. We'll write 

oo oo 

and likewise for h, then expand the nonlinear term: 



(53) 



OO OO OO 

: E (r(.g, fy), + E E« r ^ Vu), /i> + < r (-9, /j+j 



)} 



Consider the sum over I of the terms (T(g, hj+i), f 3 ) for fixed j. We expand T as a 
series by introducing the cutoff around the singularity of b in terms of and : 



oo oo oo 

E < r (-9' Vu), /j) = E E (ff . *i+ufi) - T-(g, h j+l ,fj)} 

l—l k— — oo l — l 



(54) = E {TU9,h-P j h,f j )-T'l(g,h-P j h,f j )} 

(55) 



k— — oo 



E E i T + & ^ . Si) - T - ■'■>■ h i+i >Si)}- 



1=1 k=l 



Here we have used the basic telescoping property h — Pjh = Y]?Z-f hj+i ■ Also, 
throughout the manipulation, the order of summation may be rearranged with 
impunity since the estimates we employ below will imply that the sum is absolutely 
convergent when g,h,f are all Schwartz functions. Regarding the terms (54), the 
inequalities (30) and (22) dictate that 

o 

E \T$(S, h - Pjh, fj) - T k l{g, h - P h, fj)\ < \g\ s „\h - P 3 h\nAf 3 \L* +2s - 

k— — oo 

Since \Pjh\ v ,s ;$ |^|?7,<5 ( a consequence of (51)), one may conclude that 

oo oo 

E E i T +fo h - p A Si) - t_(s, h - p^, fj)\ < e \Si\L> +2a 

j=0 k=-oo -'- n 



3=0 



< 



T},5 



uu 

E^'i/iii^ 



< 



\g\Li\h\ L 2 +2s + \g\ L * +23 \h\Lz) I/Ijv-.t. 



This is just Cauchy-Schwartz. The expansion of the product Igla^l^lr/,^ proceeds 
by (15) and the favorable comparison of the square-function norm of / to the norm 
|/|jv*.-r is provided by (52). As for the terms (55), when k < j a similar approach 
holds; namely, (31) and (24) guarantee that 

E |?t (.9, h j+l ,fj) - Tl (g, h j+l ,fj)\< 2 2s %| i2 \h j+l \ L 2 +2a \fj | 
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(we have used the trivial facts that 2 2sk = 2 2s ^2 2s( ~ k -^ and J2l=i 2 2s( - k -^ < 1). In 
particular, this inequality may be summed over j; another application of Cauchy- 
Schwartz gives that 



oo J 



J2J2\ T +(^ hj+i . h ) - Tl (g, h j+l ,fj)\ 

3=0 fe=l 



< 2~ sl 



3=0 



3=0 



< 2 



g\L*\h\N°n\f\N*n 



This estimate may clearly also be summed over I > 0. 

A completely analogous argument may be used to expand T for the terms in (53) 
of the form (T(g, hj), f j+l ) in terms of T k - T k ; 



OO OO 



1 = 1 



(56) 



(57) 



^(T^hj), f j+l )= J2( T +- T *K9,h j ,f j+ i) 

A:— — oo /— 1 


= J2 ( T +- T *)(g,h j ,f-P j f) 

OO OO 



k— — oo 



1=1 k=l 



In this case the estimates (30) and (26) are used to handle the terms (56) just as 
the corresponding terms (54) were handled. Regarding the sum (57), now (31) and 
(28) are used to estimate the sum analogously to the estimates of (55). The only 
difference is that the roles of h and / are now reversed. 

Recalling the original expansion of (T(g, h), f) it is clear that the only terms that 
remain to be considered are the following: 



OO OO OO 



(58) f( 9l hJ) d = l J2J2 £ /../•* T*(g,h j+l ,fj)} 

1=0 j=0 fc=j+l 



OO OO OO 



(59) 



+ £ £ £ i T +(9, h j} f j+ i) T k (g, hj, f j+l )} . 

1=1 j=0 k=j+l 

In other words, we have already established the inequality 
(T(g,h),f)-f(g,h,f)\ < \g\ L *\h\Nsn\f\ N s„ 

+ \9\ L * +2 „ [\h\L*\f\N:f + \h\N..-,\f\ L >]. 

Then the terms on the right-hand side of (58) and the terms (59) are both treated 
by the cancellation inequalities. The terms (58), for example, are handled by (38) 
and (39) (when s > ^). For any fixed l,j, k, we have 



\T*(g,h j+l ,fj) -T k {g,h 3+l J 3 )\ 



<2 2sk \g\ L ,\h J+l \ L , 



2=0 



ik—€ik 
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for appropriate nonnegative e-j's satisfying 2s — i — e, < (note that both (38) and 
(39) may be written in this form). In either case, there is decay of the norm as 
k — > oo since 2s — i — < 0, so 

oo 

\Tl{g,h ]+l ,f 3 ) -T k :(g,h ]+h f 3 )\ 

k=j+i 



< 2 2sj 



9\L*\h j+ i\ L z +2t 



2 

E : 

i=0 



-Jl 



v 4 |7j 



Just as before, Cauchy-Schwartz is applied to the sum over j. In this case 2^ 2s ^ 1 ^ 
is written as 2( s ~ i > 3 2 s ( J+l >2~ sl \ the first factor goes with /, the second with h, and 
the third remains for the sum over I. Once again (52) is employed. Finally, the 
factor of 2~ sl allows one to finish the sum over I. 

The desired bound for the nonlinear term is completed by performing summation 
of the terms (59). The pattern of inequalities is exactly the same as the one just 
described, this time using (45) and (46). In particular, one has that 

\Tl{g,h h f j+l ) -Tffahjjj+Ol 



<2 2sk \g\ L *\f; 



3+l\Ll 



— ik — ei k 



v 4 r/i, 



Again, both (45) and (46) may be written in this form with q > satisfying 
2s — i — £j < 0, leading to the corresponding inequality for the sum over k: 



E 

k=j+i 



Tl(g,h J+l J,) —T*(g,hj + i,fj)\ 



< 2 2s J 



9\L^\rij+i\ L 2 +23 



E 2 ~ iJ 'i v *r/, 



The same Cauchy-Schwartz estimate is used for the sum over j\ there is exponential 
decay allowing the sum over I to be estimated. The end result is precisely (16). 

5.2. Upper bounds related to the linear operator. In this section, we prove 
Lemma 3 and Lemma 5. The inequality (17) follows by estimating (T(M, /), /) us- 
ing (16) (plus the basic observation that \/\l 2 ^ \ f\N a >~>)- Regarding the "compact" 
piece of the linear term (i.e., the operator (14) and associated inequality (18)), one 
would like to establish the inequality in Lemma 5. 

To achieve Lemma 5, we essentially reprise the arguments of the previous section, 
in a slightly simpler form. Now M coincides with e~ V4 ^ 2 restricted to the paraboloid 
(v, ^M 2 ) as has been already noted. In particular, then 

M~ 1 P j e l (v)= t dv' {v')2 3j (p{2 j {v-y!_))p{y!_), 

where <p(v) = e 2 3 lvi (p(v). Recall ei(v) is defined above Lemma 5. Thus, the 
asymptotic expansion (49) applied to (p yields that 

^{ v + 2 3 t v u - v)) = ip(z v u) + 2"V (1) (r„u) + Ej(u, v), 
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for some odd function ip 1 -- 1 ^ . Thus, just as was done for (50) (namely, expanding (p, 
p and {%)') and noting that all of the terms of order 1 and 2~' J vanish by symmetry 
leaving only terms of order 2 _2 - J and higher), we may deduce that 

\M- 1 Q j e l {v)\<2-V (v) A , 

and, more generally, 

|A/- 1 |V4| m Q J e ; W| <2(- 2 +™)» 4 , 

for to = 0, 1, 2 (note that the power 4 arises from the fact that (v) < (v) 2 and the 
polynomial is quadratic). Informally, these inequalities guarantee the intuitively 
obvious point that the function ei is, in fact, very smooth as measured by our 
Littlewood-Paley projections (in fact, the decay rate 2 ly ~ 2+m ^ is limited only by 
the fact that we have not required that the scaling function ip have many vanishing 
moments, in which case the decay would be even better) . Now we expand the sum 

oo oo 

(T(g, e,), })= EE i T +(9, ei h f) - T*{g, ey, /)} , 

k= — oo j=0 

where we abuse the future notation by defining e;o == Po e i and = f Qjei for j > 0. 
To estimate the summand, we use both (34) and (27) on and T* separately when 
k < j, and when k > j we use the cancellation inequality (48). In both cases we also 
use the estimate for AI~ 1 Qjei to obtain an inequality for 1 1 V 1 e / j I i 2 + 2 i namely, 

that it is bounded by a uniform constant times 2^ _2+m ^ . Thus our inequality for 
the sum becomes (after exploiting (34), (27) and (48)) 

oo oo 2 

l(r(ff,e,),/>| < \9\ v> 6\fkv E E E min{2 2fcs 2-^,2( 2s - m ) fc 2(- 2+m W2- e - fc }, 

fe=— oo j=0 m=0 

for appropriate nonnegative e m 's which satisfy 2s — to — e m < 0. Clearly this 
sum must be finite, which establishes the desired inequality for the compact piece. 
Moreover, since each estimate for a trilinear form of (g, h, f) has a corresponding 
analog with the roles of h and / reversed, an identical argument to the one just 
given (using the estimates (30), (22) and (47)) establishes that 

\(ngJ),ei)\<\gWs\fkv, 

holds uniformly for r\ and S as well. In particular, Lemma 5 follows. 

6. The main coercive inequality 

This section is devoted to the proof of Lemma 4. Our approach involves direct 
pointwisc estimates of a Carlcman representation in Section 6.1. However this 
argument will not be completely sufficient, as explained below. Thus in Section 6.2 
we prove an estimate dubbed "Fourier redistribution" to finish the desired bound. 
Then in Section 6.3 we will establish functional analytic results on the space iV s ' 7 . 

6.1. Pointwise estimates. For any Schwartz function /, consider the quadratic 
difference expression arising in the study of (Nf,f) from (13). By virtue of the 
Carleman-type change of variables, it is possible to express this semi-norm as (20), 
where the kernel can be computed with Proposition 11 in the Appendix to be 

(KM Iff M def o/ A M * M * o fo I ' V '-< \ 

(60) K(v,v)=2 dm v > — -r B [2v - v - v„, -— . 
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Above M* = M(v'^ + v' — v) and the integration domain is the hyperplanc 

^=Ket 3 I (t/ -«,<-«) =o}. 

Then denotes the Lebesgue measure on this hyperplanc. 

Our goal is to estimate this kernel K pointwise from below and compare it to 
the corresponding kernel for the norm | ■ from (8); this, by virtue of (52), 

allows control of our anisotropic Littlcwood-Paley square function by (Nf, f) . We 
make this estimate when \v — v'\ < 1 and \\v\ 2 — \v'\ 2 \ < \v — v'\. This constraint 
will require the introduction of a somewhat technical argument, but it is necessary 
since the required pointwise bound fails to hold uniformly outside this region. 

On the hyperplanc E%,, we have (u — v', V — =0 and \2v — v'—v'J\ = \v' — i>*|; 
in particular, then 

2v — v' — vt v' — v 



[2d — v' — vl\ \v' — v 
By virtue of the lower bound for &(cos#) in (4), it follows that 
/ v' - v' \ W - v' l 2+2s 

B ( 2 - - <. y^n) Z *(!«' - <D |,_,;|U V<i>k4 

The indicator function must be included because of the support condition in (4). 
Thus the kernel K(v, v') from (60) is bounded below by 

(61) | V - v'\- 3 ~ 2s I dir K MM'MW ~ <\)W ~ <f +2s l|,-<|>M|- 

Next we consider the magnitude of the projections of u* = v' + v' t — V and v'^ in the 
direction of v — v' . The orthogonality constraint (v — v' ,v — v' t ) =0 dictates that 

v — v'\ I i v — v'\ — \v — v'\ 2 + \v\ 2 — \v'\ 2 



\v — v'\ / \ ' \v — v'\ / 2\v — v' 

v — v'\ I v — v' \ \v — v'\ 2 + \v\ 2 

V, ■ 



\v — v'\f \ \v — v'\ I 1\v — v'\ 

With our assumptions \v — v'\ < 1 and |M 2 — \v'\ 2 \ < \v — v'\, both right-hand sides 
are uniformly bounded by 1 in magnitude, implying that \v*\ 2 + |t^| 2 < 2|u^| 2 + 1, 
where w'^ is the orthogonal projection of v'^ onto the hyperplanc through the origin 

with normal v - v' , e.g. = v' t - j^pj (^F]> v '*\ 

This implies M*M* > e - ' 1 "*' I 2 uniformly. Let w' and w be the orthogonal 
projections of v' and v respectively onto this same hyperplane through the origin 
with normal v — v' . Trivially \v' — v'^\ > \w' — w'J\. Further \w' — v'\ < 1 since 

\(v',v-v')\ = i \\v\ 2 - |u'| 2 -|u-u'| 2 | < \v-v'\. 

Write v'x — + v — w, then we may parametrize the integral in (61) as an integral 
over uu'^ (with unit Jacobian) and thereby bound (61) uniformly from below as 



K(v,v') > \v-v'\- 3 - 2s I d7r. <e -^*lV ~w'^ +2s+1 

E' 



with E' = {w'x | (w'^jV — v') — 0, \w' — w'J > \v — v'\}. This change of variable 
preserves these estimates because the effect of parametrizing the integral with w' t 
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is to shift the plane in the direction v — v', while the inequalities for M*M^ and 
\v' — v^l are in terms of only the components orthogonal to v — v' . 
If | w' | < 4, then since \v — v'\ < 1 it is not hard to see that 

dir w , e-^KlV - w'^ +2s+1 > f dir w , e~^<f > 1. 

JE'n{\w'-w',\>i} 

When \w'\ > 4, we may restrict w'^ to lie in the disk ^\w'\ > |iy^| + 1, which implies 
in particular \w' — w' t \ ~ \w'\. We thus have the following: 

r h\w'\-l 

dir w ,e-iW 2 \w' - w'^ +2s+1 > ( w ') 7+2s+1 / dp pe~^ 2 > ( w ') 7+2s+1 . 

E' JO 

Since \w' — v'\ < 1, the final, uniform estimate for (60) becomes: 

KM > \v-v'r 3 ~ 2s (t;')' y+2 * +1 i| 3 ,- fi i|< 1 l||H»-i«'i a l<l»-« 'I- 

On this region \v — %/\ < \v — v'\ and (v) w (t/), so with (13) we have uniformly 

(62) |/||>/ dv f dv> {f ~^ 2s ((v) (i/» 3±¥±1 l| t ,-«,/|<il|| t ,|»-y|»|<|w|- 

To obtain a favorable coercivity estimate from this, it would suffice to show that the 
expression (62) is bounded from below by the corresponding piece of (8) (since the 
former expression has already been shown to be connected to our exotic Littlewood- 
Paley projections). Because of the cutoff restricting l|M2_|.u/| 2 |<b— u'l a direct point- 
wise comparison is not sufficient. This is not merely a limitation of the argument 
leading to (62); in fact, a more involved analysis of (13) shows that there is ex- 
ponential decay of K(v,v') in \v — v'\ when v and v' point in the same direction. 
Thus there is an intrinsic obstruction to obtaining the correct coercive inequality 
by means of a simple, pointwise comparison of these expressions. 

6.2. Fourier redistribution. To get around this obstruction, we use the following 
trick (dubbed here "Fourier redistribution"). Essentially the idea is to appeal to 
the Fourier transform in the situation where the pointwise bound is not available. 
The key idea is already contained in the following proposition: 

Proposition 8. Suppose K\ and Ki are even, nonnegative, measurable functions 
on R 3 satisfying 

|2 



du Ki{u)\u\ < co, i = l,2. 

Suppose 4> is any smooth, nonnegative function on R 3 and that there is some con- 
stant C<f, such that |V 2 </>(u)| < for all u. For I = 1,2, consider the following 
quadratic forms (defined for arbitrary real-valued Schwartz functions f): 

I/Ik-, = I dv f dv' <j>(v)W)Kx(v v%f(v) f{v')) 2 . 

If there exists a finite, nonnegative constant C such that, for all (€l 3 

/ du K^u)^ 1 ^ - 1| 2 < C + f du K 2 (u)\e 2 ^' u) - 1| 2 , 

then for all Schwartz functions f, 

\f\ 2 Kl <\f\ 2 K2 +C'C^ [ dvcf>(v)(f(v)) 2 , 
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where the constant C satisfies C < 1 + C + J R3 du{K\{u) + K2{u))\u\ 2 uniformly 
in K\,K2, 4> and C. 

Proof. We begin with the following identity: 

<f>{v)<f>{v'){f{v) - f{v')f = (^v)f(v) - 4>{v')f{v')f 

+ ^v)(f(v)) 2 (^v')- ( j ) (v)) 

+ <l>(v , )(f(v')) 2 (<f>(v)-cf>(v')). 

Multiply both sides by Ki(v — v') and integrate with respect to v and v'. Exploiting 
symmetry, the result is: 

dv f dv'K^v - v')(f(v) - f{v')fcf>(v)ct>{v') 

R3 JR 3 



dv I dv'Ki(v - v')(</>(v)f(v) - ^{v')f{v')) 2 

3 

dv(j)(v)(f(v)) 2 p.v.f dv' Ki{v - v'){4>{v') -4>{v)). 



Now Taylor's theorem and the hypotheses on the second derivative of <f> dictate 

W) - 0(«) - («' - v, V<Kv)) \<\\v'- v\ 2 C^. 
If we define C(Ki) = / duKi(u)\u\ 2 , it follows that the difference 
dv [ dv'K^v - «')(/(«) - f(v')) 2 <P(v)<f>(v') 

dv f dv'K^v - v'){<P{v)f{v) - </>(v')f(v')) 2 



is bounded above by C(Ki) J R3 dv<p(v){f (v)) 2 . If we cutoff \u\ > e, then clearly the 
Plancherel formula can be applied to the second term inside the absolute values 
above, with F(v) = <p(v)f(v), to get 

dv [ du Ki(u)(F(v +u) - F{v)) 2 l H>e 

JR 3 

(63) = / d£ [ duK^e 2 *^ -l\ 2 \F(0\ 2 l H>i . 

JR 3 JR 3 

Clearly the limiting case e — > will hold as well because \e 2 ^ 1 ^^ — 1| 2 vanishes 
to second order in u and F may be assumed to have arbitrarily rapid decay in |£|. 
From here, the remainder is clear. The hypotheses on K\ and K2 give that the limit 
of the Plancherel term (63) is bounded above by the Plancherel term for K2 plus 
C times the L 2 -norm of F. This term plus the errors in comparing the Plancherel 
pieces (63) to the norms | • \ 2 Ki give rise to the constant C. □ 

Next, fix functions Ki, Ki on E 4 given by K\(u) = 3 ~ 2s l| li |< 1 and K2{u) == 
|u| _3_2s l| u i <1 l| U4 | <e | M | ! that is, K\ is restricted to the unit ball and K2 is further 
restricted to (1 — e 2 )u 2 < e 2 (u\ + u\ + u 2 ). We define the semi-norm Nq by 

l/lk=/ dv I dv'K 2 (v-vO(f'-f) 2 ({v'){v)r 
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Note that, if K 2 is replaced by K\, the resulting expression is the derivative part 
of our main norm (8). By a pointwise comparison of K\ and K 2 , it is trivially true 
that |/|jv I/I-/V 3 'T, but our goal is to prove an inequality in the reverse direction. 
To that end, let {4>} be a smooth partition of unity on R 4 which is locally finite 
and satisfies uniform bounds for each <p and their first and second (Euclidean) 
derivatives. Suppose furthermore that each <fi is supported on a (Euclidean) ball of 
radius |. Recall the notation from Section 4.1. We restrict these functions to the 
paraboloid (v, ^\v\ 2 ) and insert them into the norms | • \n',i and | • |jv : 

(64) [ dv f dv'K l (v- 2 0(f-f) 2 (v) : ^(v') : ^<l>(y)<t>{£), 

for I = 1,2. Suppose that vq £ M. 3 satisfies 4>(vq) 7^ for some fixed </>. Make the 
change of variables uh>«o + t Vq u and likewise for v'; including the Jacobian factor 
(vq) for each integral, the result is an integral over u and u' of the integrand 

(v y 2 Ki^VqJtJ^U - Vq + T Vo u' )(f(v a + T Vg u) - f (l) + T Vo u')) 2 

x <f>( vo + t Vo u )4>( v + T Vo u' )((v + t Vo u) (v + t Uo m'))" + 2 + . 
Now we expand. The argument of Ki, for example, becomes 

Vq + T Vo U - Vq + T Vo u' = T V(j (u - u') + ^{\u\ 2 - K| 2 )&4 

, A / . u + u'\ 
= Z vo {u-u) + (u-u, ) e 4 . 

Since the support of (f> is in a ball of radius | , it follows that | \ < | , hence the 
magnitude of the coefficient of is at most f |u— u'\ = | \r v (u — u')\, so 

\( VQ + TypU - VQ + T va u! ) -T V0 {U-U )\ < ^\t Vq (u — u')\. 

Thus on this piece of the partition, K\ may be bounded above by 

K^ Vp + T Vg U - V + T Vo u' )4>( v + T Vo u )<j){ vo + T VQ u' ) 

< \u - u'\^ 3 ~ 2s l lu _ u , l < 2 4>( v + T Vo u )(j)( vo + T Vo u' ), 

which is translation- invariant in u and u' . 

Next we must make a similar bound for K2 from below. On this piece of the 
partition, if the fourth coordinate of t v (u ~ u') is bounded above by | |t„ (u — it')|, 
then the fourth coordinate of (vq + t Vq u — vq + t Vq v!) will be at most 



^\T V0 (U - U')\ + ^\t Vo (u - U')\ < - \ Vq + T Vo U - Vq + T VqU ' \. 

Thus it also holds that 

\u - ^r 3 ~ 2s l| M -„'|<i l(„ o )-i|(^ o , M -„')|<j|„-^|0( «o + T Vo u )<j>( v + t Vo v! ) 

< K 2 ( vq + T VQ U - Vq + T Vo u' )(j)( v + T Vo u )(j)( v + T Vo u' ), 

which is also translation-invariant. To apply the proposition, then, it suffices to 
check the Fourier condition and estimate the derivatives of the cutoff functions. 
Clearly zeroth-order through second-order derivatives of 

-1 + 23 + 1 



(vq) (vq + U) 2 (j){V(j +T V0 U), 
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-1 + 23-1 

with respect to u will be uniformly bounded by (vq) 2 by virtue of the corre- 
sponding estimates for <f> coupled with the fact that t Vo has norm 1 as a mapping 
of Euclidean vector spaces and r Vg is an isometry. 

Modulo the verification of the Fourier condition, then, the norm |/||/- s , 7 on any 
particular piece of the partition will be bounded above by \f\% on the same piece 
plus an error term which is given by integration J du(f(vo + t Vq u)) 2 4> times a con- 
stant linear in <fi. That is, the difference of the quantities given by (64) for I = 1 
and / = 2 is at most 

C'Ct / du cf)(v + T Vo u)\f(v a + t Vo u)\ 2 , 

as a result of the previous proposition. Here <j) = f (■y)( 7+2s ~ 1 )/ 2 ^ (recall that the 
extra factor of (vo) comes from the change-of- variables we employed). Thus the 
quadratic dependence on </> gives a factor of (vq) to the power 7 + 2s — 1; however an 
additional factor of (vo) is obtained when the change-of-variables is reversed (that 
is, Vq + T V( ,u reverts back to v). Thus, summing over the partition will give 

\f\ 2 N ^<\f\k+ [ dv(f(v)f ( V y +2s . 

JR 3 

To complete the comparison, then, it suffices to make the following estimate: 

Proposition 9. Fix any e > 0, and let E\ and E2 be the sets in M 3 given by 
E x = {ueR 3 I \u\ < 2 } and E 2 = {u G K 3 | |«| < \ and |u 3 | < e\u\ }. Then 



(65) / du | e 2 ™<«>"> - l| 2 | u |- 3 - 2s < 1 + / du | e 27 "<«-«> - i n 

J Ei 

uniformly for all (el 3 . 



2|„,|-3-2s 

j, a \c •- ■ — ±| 

E2 



Proof. Writing both sides in polar coordinates, we see that each side may be realized 
as an integral over the unit sphere S 2 of 

J Ei 

where E\ = § 2 , E2 is a small band near the equator, and 'I' (A) is of the form 

*(A) = f f dt |e 2wiAt - II 2 *- 1 " 2 *, 
Jo 

for some appropriate value of a (a = 2 or a = ^). From the elementary inequalities 



L 



2 • 

(2A)- 1 A2X)- 1 

dt \ e 2*i\t _ i^-l-a, „ / dt> 2 t 2 t -l-2s ^ X 2s 

JO 

/■oo 

2-wiXt t 1 2 . — 1 — 2s I j. . — 1— 2s ^ \2s 



dt\e lmxt -\\ z t-'- ls < \ dtt- L ~ As ^X 2 

(2A)- 1 J(2\)~ 1 

it follows that the integrands will be comparable to | (£,<r) \ 2s when this quantity 
is bounded below by a fixed constant and less than a constant times | (£,cr) | 2s 
regardless of whether or not this quantity is bounded below. For any £ with |£| > 1, 
then, at least a positive measure region of § 2 will have | (£, a) \ > |(| (whether in 
Ei or E2), so both sides of (65) will be comparable to |£| 2s , which is sufficient for 
the inequality (65) to hold. □ 
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The proof of the coercive inequality is now complete, for we demonstrated that 
dv f dv, I da B(f - f) 2 MiM*>\f\ 2 NQ , 

JR3 JS 2 

by direct pointwisc comparison and that \f\% + \f\ 2 L 2 <^ |/|w s ^ by Fourier 
redistribution. The combination of these inequalities gives Lemma 4. 

6.3. Regarding the functional analysis of N sn . An important consequence of 
the analysis of the previous section is that we have an alternate characterization 
of the space N sn in terms of the usual Sobolev spaces. In particular, let {</>.;} 
be a partition of unity constructed as above by restricting a smooth, locally finite 
partition of unity on K 4 (such that each (pi has support in a ball of unit radius) to 
the paraboloid (v, ^\v\ 2 ). For each 4>i in the partition, let Vi be some point in its 
support. If wc define 

fi(u) = f (f>i( Vj + T Vi u )f(vj + r Vi u), 
it follows that we have the comparison 

oo 

(66) I/Ia^«E^> 7+2S_1 |/^' 

i=l 

where H s is the usual (three-dimensional) L 2 (R 3 )-Sobolcv space. This result is 
true by virtue of the fact that 

1/1*. ~ \^ + JJ V jj V | w _ w /|3 +2 . Wl<X. 

which follows itself by an application of the Plancherel theorem as in Proposition 
8 together with the asymptotic estimates for the integrals (65). 

With the aid of (66), a number of elementary functional analysis properties of 
N sn reduce to the situation of the standard Sobolev spaces. For example, it is a 
simple exercise to show that Schwartz functions are dense in 7V S ' 7 by exploiting this 
same fact for the space H s , approximating individually in H s , and summing over 
the partition (note that this requires the elements of the partition <j>i themselves to 
be Schwartz functions, but this additional restriction is not a problem to satisfy). 

A somewhat more sophisticated result which may be obtained by similar rea- 
soning is the fact that N sr< embeds compactly in L 2 . Clearly one also has the 
comparability statement that 



E (^r 1 i-ftii 2 = X 3 dv (e i-m^i 2 ) i/mi 2 , 

and the sum of the squares of the <?Vs must be uniformly bounded above and below. 
Now suppose /" is any sequence of functions in N s ' 7 with norms uniformly bounded 
by 1, having the weak limit f° E TV''' 7 . Then the compact embedding of H s into 
L 2 on compact domains ensures that /" converges in L 2 for each v. the limit of /" 
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:»,<) 



must equal ff. Now by (66) and the corresponding estimate for the L 2 -norm, we 
have for any R > 1 that 

oo 

ir-zV^E^r 1 !^-/?!^ 

i=l 

< E ^r 1 m ~ + r-"- 2s E ^) 7+2s_1 iff - f?\n- 

i ■ \vi\<R i : \vi\>R 

< E K)" 1 i/r-/rii=+^" 7 " 2s ir-/V— 

i ■ \vi\<R 

To show that f n converges to f° in L 2 , simply observe that the finite sum over i in 
the last line above goes to zero by compactness for any fixed R, while \ f n — f°\N s --i ^ 
1 uniformly for all n; thus taking R — > oo establishes the claim. 

7. DE-COUPLED SPACE- TIME ESTIMATES AND GLOBAL EXISTENCE 

In this last section, we show that the sharp estimates proven in the previous 
sections can be applied to the modern technology from the linearized cut-off Boltz- 
mann theory to establish global existence. This works precisely because of the 
specific structure of the interactions between the velocity variables and the space- 
time variables. The methodology that we employ essentially de-couples the required 
space-time estimates that are needed from the new fractional and anisotropic de- 
rivative estimates which are shown in the previous sections. 

The method that we choose to utilize in this section goes back to Guo [32]. A key 
point of this approach is to derive a system of space-time "macroscopic equations," 
see (79) - (83) below, which have certain elliptic structure and also some hyperbolic 
structures. These structures can be used to prove an instantaneous coercive lower 
bound for the linear operator L, for solutions to the full non-linear equation (9), 
in our new precise norm (8). This original method [32] made use of high order 
temporal derivatives, which we could also utilize. But as a result of advances 
in [33], [23,35] the need for temporal derivatives was removed from the method. 
The key point here is to use both the macroscopic equations (79) - (83) and the 
conservation laws (84) - (86) to remove the need to estimate time derivatives with 
an "interaction functional" that is comparable to the energy. We point the readers 
attention to the general abstract framework of [52] also in this direction. 

In what follows, we will show that the use of our new precise weighted geometric 
fractional derivative norm (8) and our crucial sharp estimates stated in Lemma 2, 
Lemma 4, and Lemma 5 can be combined with the above general de-coupled energy 
method in order to prove global existence and decay as in our main Theorem 1. 

We will now discuss the coercivity of the linearized collision operator, L, away 
from its null space. More generally, from the H-thcorem L is non-negative and for 
every fixed (t, x) the null space of L is given by the five dimensional space 

(67) Kf = span{y/fi,v 1 y/Ji,V2y/Ji t V3y/Ji, \v\ 2 ^/Jl} . 

We define the orthogonal projection from L 2 (Rf,) onto the null space AT by P. 
Further expand P/i as a linear combination of the basis in (67): 

(68) Ph= | a h (t,x) +J2tf(t,x) Vj +c h (t,x)\v\ 2 1 JJi. 
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We can then decompose f(t, x, v) as 

/ = P/ + {I-P}/. 

We characterize in Lemma 6 the functional properties of the linearized collision 
operator that are useful for our main results. 

Lemma 6. L > 0. Lh = if and only if h = Ph. And 36q > such that 

{Lh,h)>S \{I--p}h\ 2 Nsn . 

The proof will follow directly from Lemma 5 and the splitting L = N + K. 
The coercive estimate in Lemma 6 is proven via a contradiction argument. This 
is the only location in our paper where a non-constructive argument is used. We 
expect that our estimates can aid in a future constructivc-but perhaps substantially 
longer proof of this coercivity; see [41,42] and the references therein. 

Throughout the remainder of the paper, we will abuse notation and write |/i|^-»,t 
from (8) in place of (Nh,h) from (12); since these quantities have already been 
shown to be comparable, all of the previous inequalities involving | ■ \n b ^ will re- 
main true after this redefinition at the price of a fixed multiplicative constant. 

Proof of Lemma 6. Most of this lemma is standard, see e.g. [27] for proofs of the 
first statements in the cut-off case. Without cut-off the first statements can be 
established with the same proofs as in the cut-off situation. This is done via the 
usual approximations of the singular kernel (4) with a non-singular kernel b e (cos 9) 
and sending e 1 0- We only prove the coercive lower bound for the linear operator. 

Assuming that coercivity fails grants a sequence of functions h n which satisfy 
Ph n = 0, \h n \ 2 Nsn = (h n ,h n ) N s,-, = (Nh n ,h n ) = 1 and 

(Lh n ,h n ) = \h n \%s n - (Kh n ,h n ) < -. 

71 

Thus {h n } is weakly compact in 7V S ' 7 with limit point h°. By weak lower-semi 
continuity \h |jy-*,-y < 1. Furthermore, 

(Lh n ,h n ) = 1 - (Kh n ,h n ). 

We claim that 

lim (Kh n ,h n ) = (Kh",h a ). 

n—>oo 

The claim will follow from the prior Lemma 5. The claim implies 

= 1- (Kh a ,h°). 

Or equivalcntly 

(Lh°,h ) = \h°\ 2 N ^-l. 

Since L > 0, we have = 1 which implies h° = Ph°. On the other hand 

since h n = {I — P}/i n the weak convergence implies h° = {I — P}/i°. This is a 
contradiction to |/i |^ s ,-, = 1. 

It remains to establish the claim. We expand out 

(Kh n , h n ) - (Kh°, h°) =(Kh°, (h n - h")) 

+ {K(h n -h°),h°) 

+ (K{h n -h°),(h n -h )). 
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We recall the definition (14) and the general estimate in Lemma 5; setting 77 = 6 = 1, 
it is easy to see that 

\(Kg,h)\<\g\ hl \h\ 1A . 

Since \f\ v ,s ;$ I/Iat^ for any fixed 77, 5, it follows that K is bounded on N sn . In 
particular then 

(Kh°, (h n - h )) -> 0, 

as n — > 00 by weak convergence. Since (A' g, h) = (fj, if h), the same is true of the 
second term. Finally, by (18), we have uniformly for all 77 > that 

I (K(h n - h°), (h n - h )) I < r,\h n - ^II^crs) + C(r,)\h n - hP\% . 

Since \h n — h°\ L 2^ < \h n — /i°|jv».t < 2 for all n and N sn embeds compactly in 
L 2 V as in Section 6.3, it follows that for any 77 > 

\imsup\(K(h n - h°),(h n - h°))\ <?7. 

n— >oo 

In particular, this implies the limit is zero and establishes the claim. □ 



7.1. Local Existence. Local existence has been shown in the recent preprint [8], 
with higher regularity assumptions on the initial data than we consider. We estab- 
lish the a priori bounds for local existence herein using our estimates above with 
initial data in the space L^H^ . Our local existence proof for (9) is based on a 
uniform energy estimate for an iterated sequence of approximate solutions. 

Our iteration starts at f°(t, x, v) = f fo(x, v). We solve for f n+1 (t, x, v) such that 

(69) (d t +vV x + N)r+ 1 +Kf n = r(f n r+ 1 ), r +1 (0,x,v) = Mx,v). 
For notational convenience during the proof we define the "dissipation rate" as 

*>(/(*))= E ii 5a /(*)iiW 

|a|<7V 

We will also use the following total norm 

(70) £?(/(*)) = ll/(*)lli'(R3 ; ^(T3)) + J dr 2?(/(r)). 

Our goal will be to obtain a uniform estimate for the iteration on a small time 
interval. The crucial energy estimate is as follows: 

Lemma 7. The sequence {f n (t,x,v)} is well-defined. There exists a short time 
T* = T*{\\fo\\ 2 L2HN ) > 0, such that for \\fo\\ 2 L2H N sufficiently small, there is a 
uniform constant Cq > such that 

(71) sup sup g(f n (r))<2C \\fo\\ 2 L 2 H N. 

n>0 0<t<T' v x 

Proof. The proof proceeds with an induction over k. Clearly k = is true and we 
assume that (71) is valid for k = n. For a given f n , it is standard to show that 
there exists a solution f n+1 to the linear equation (69). We focus here on the proof 
of (71). Take the spatial derivatives d a of (69) to obtain 

(72) (ft + v ■ V,) d a f n+1 + N (d a f n+1 ) + K (d a f n ) = d a T (/", f n+1 ) . 
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Therefore, applying the non-linear estimate in Lemma 2 yields 

^ii^r +1 iii;^ + ii« a r +1 iiij JV ^ + (*'(« a r) ) s a / n+1 ) 

= (d a T(f n f l+1 ) ,d a f n+1 ) 

~ L?JI/" + 1 ||fl-«7V=.^ + ll/ n |lffi v A rs '-'ll/ ,l+1 |lff^iV s .T||/ n+1 ||H«L2- 

Then integrating the above over [0, t] we obtain 

hdV n+1 (t)\\l 2L2 + f dT\\d a f n+1 (r)\\l 2NS , 7 + f dr (K(d a f n ),d a f n+1 ) 
1 " x Jo Jo 

(73) < hd a fo\\h vLl +C f dr \\n H »L>\\f n+1 \\ 2 H»NsAT) 

+C f dr \\r\\H-N^\\f n+1 \\H-NsM\f n+1 \\H-L^r). 
Jo 

We notice that from Lemma 5 applied to (14), for any 77 > small and 5 = 1/2, 
dT(K(d a f n ),d a f n+1 ) 



1 '1, 



< / dr[-\\dV n+ \T)\\l 2+2s +CJdy n+ \r)\\l l 



+ V dr\\d a r(r)\\i 2 +C dr\\d a r{r)\\h ■ 
Jo " l+2s Jo 

We incorporate this inequality into (73) and sum over |a| < N to obtain 

g(r +1 (t))<Co\\fo\\i lH K + [ ^ {c\\r+Y H »Li(T)+cv\\r\\ 2 H »Li +2 >)} 

J 

+0, / rfr + C sup S(/" +1 (r)) sup e 1/2 (/ n (r)) 



0<T<t 0<T<t 



< Co || /oil + C„^ sup 5(/ n+i (r))+ sup g(/ n (r)) 

" * L0<r<i 0<T<t 

+C V sup 5(/ n (r)) + C sup Q{f n+1 {r)) sup ^(/"(r)). 

0<r<t 0<r<t 0<T<t 



We are using the total norm from (70). By the induction hypothesis (71) 

2 



sup G(f n (T)) < 2Co||/o|| 2 



0<T<t 

Then we collect terms in the previous inequality to obtain 

{l-<^T*-a||/o|U ?i ^} sup G(f n+1 (t)) 

0<t<T* 

< {Co + 2C V + 2C V T*C } \\fo\\ 2 L 2 H N. 
By choosing rj small, then choosing T* = T*(\\fo\\ L 2 H N) small, we have 

sup G(f n+1 (t)) <2C \\fo\\h H ,. 

0<t<T* " x 

We therefore conclude Lemma 7 if T* and ||/o|| 2 2#jv are sufficiently small. □ 

With our uniform control over the iteration from (69) proved in Lemma 7, we 
can now prove local existence in the following Theorem. 
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Theorem 8. (Local Existence) For any sufficiently small Mq > 0, there exists a 
time T* = T*(M ) > and Mi > 0, such that if 

\\M\l lH » < Ah, 

then there is a unique solution f(t,x,v) to (9) on [0, T*) xT^xR^ such that 

sup g(/(t))<M . 

0<t<T* 

Furthermore G(f(t)) is continuous over [0,T*). Lastly, we have positivity in the 
sense that if Fo(x, v) = fj, + [J, 1 ' 2 fa > 0, then F(t, x, v) = fj, + ^l 2 f(t, x, v) > 0. 

Proof. By taking n — > oo, we have shown sufficient compactness from Lemma 7 to 
obtain a strong solution f(t, x, v) to the Boltzmann equation (9) locally in time. 

To prove the uniqueness, we suppose that there exists another solution g with the 
same initial data satisfying sup < r <y» Q(g(r)) < Mq. The difference f — g satisfies 

(74) {d t + v • VJ (/ — g) + L (f -g)=T(f- g, /) + L (g, f - g) . 

We apply Lemma 2 and the Sobolev embedding H 2 (T^.) C L°°(T^.) to obtain 

l({r (/ -g, /) + r ( g , f- g )},f- g )\ < {hh^ + \\fhi m } 11/ - g\\%..-, 

+ {\\9\\h*N^-' + \\f\\H*N°--'} \\f - g\\N°.v\\f - g\\L*_ x - 
The Cauchy-Schwartz inequality (applied in the time variable) shows us that 

dr {\\g\\ H gN»>-y + \\f\\HiN°.i} \\f - g\\N°,y\\f - ghi^) 

1/2 



< VM sup ||/(r) - 5 (r)|j£ 2 / dr ||/(r) - g{r)\\ 2 Nsn 

\Q<T<t • Jo 

We have just used the following fact, which follows from the local existence, that 



Sup \\f(T)\\ L 2 H 2 + / dr ||/(t)||^ 2JV3 , 7 < M . 

0<T<t Jo 

And similarly for g(t). Since L = N + K, we use Lemma 5 applied to (14) to obtain 

(L{f - g), f-g)> l -\\f- gf Nan -C\\f- gf Llx . 
We multiply (74) with f — g and integrate over [0, i] x Tj x Rl to achieve 

\\W) 9(t)\\k,„ + \ f dr ||/(r) - g(r)\\%.„ 

2 



<y/M sup ||/(r)- 5 (r)||i a +/ dr ||/(r) - <?(r)||^ 

\0<T<t 



+ / dr ||/(r)-ff(r)||ij . 
Jo 

We deduce f = g and the uniqueness from the Gronwall inequality. 
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To show continuity of Q(f(t)) in time, we sum (73) over \a\ < N and integrate 
from ti to ti (rather than over [0, f]). Then with /" = f n+1 = f we obtain 



\g(.f(h))-g(f(h))\ = 



\\\m)\\l iH » - l\\f(h)\\l lH N + dr 2?(/(t)) 



<{l+ sup VGTmj} f 1 dr ||/(r)||W,, 

L t 2 <r<ti J Jt 2 

as ti — > t% since ||/(i")||^Arj Vs , 7 is integrable in time. 

Wc now explain the proof of positivity. Previous works which obtain the posi- 
tivity of solutions without cut-off, to our knowledge, are only [49] and [8]. We use 
the argument from [8], however their initial data is somewhat smoother than ours, 
e.g. they effectively work in /o £ H£? v for M > 5, since Fq = fi + ^/Jlfo, and they 
study moderate angular singularities 0<s<l/2. If our initial data is in H^f v , 
then since we have proven the uniqueness, we conclude that F = /i + J]If > if 
initially Fq = fi + ^/Jlfo > 0. The argument is finished by using the density of H% v 
in the larger space L 2 V H^ (T^, x K.^), standard approximation arguments, and our 
uniqueness theorem. For the high singularities, 1/2 < s < 1, the positivity can be 
established by using high derivative estimates fit) € H^f v from [29], and following 
the same procedure [8] as in the low singularity case. □ 

7.2. Coercivity estimates for solutions to the Non-Linear equation. The 

following is a by now well known statement of the Linearized H-Theorem [32] ; wc 
prove it for the first time in the regime where there is no angular cut-off, e.g. (4). 

Theorem 9. Given the initial data f € L 2 : H N (T|)) for some N > 3, which 
satisfies (7) initially and the assumptions of Theorem 8. Consider the corresponding 
solution, f(t,x,v), to (9) which continues to satisfy (7). 
There exists a small constant Mq > such that if 

(75) ||/(t)|| \ lHS < M , 

then, further, there are universal constants S > and C2 > such that 
J2 \\{l-P}d a f\\ 2 NS Jt)>5 J2 \\Pd a f\\ 2 NS Jt)-C 2 ^, 

\a\<N \a\<N 

where X(t) is the "interaction functional" defined precisely in (93) below. 

Wc prove this theorem by an analysis of the macroscopic equations and also 
the local conservation laws. The system of macroscopic equations comes from first 
expressing the hydrodynamic part P/ through the microscopic part {I — P}/, up 
to the higher order term T(/, /) as 

(76) {d t + v v x }Pf = -d t {i - p}/ + /({i - p}/) + r(/, /), 

where 

(77) /({I - P}f) d = f -{v ■ V x + L}{I - P}/. 



Notice that we have isolated the time derivative of the microscopic part. 
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To derive the macroscopic equations for P/'s coefficients a^(t,x), b{(t,x) and 
cf(t,x), we use (68) to expand entries of left hand side of (76) as 

3 

^ {vidic\v\ 2 + {d t c + dib t }vf + {d t h + d l a}v i } ^[Ji 
i=i 

3 

+ ^ y^{djbj + djbijviVj y/JI + d t a ^ffi, 

i=i j>i 

where di — d Xi above. For fixed (t,x), this is an expansion of the left hand side of 
(76) with respect to the following basis, {ek]]? = i, which consists of 

(78) (",H 2 ^)kk 3I K 2 v^)i<,<3' (^jVM)i<i<j<3 . KVm)i<*< 3 < v 7 ^- 
From here one obtains the so-called macroscopic equations 



(79) V x c = -d t r c + l c + T c 

(80) dtc + dibi = -dtn + k + Ti 

(81) dibj + djh = -dtrij+kj+Tij (i ^ j) 

(82) dth + d t a = ~d t r bl + l u + T bl 

(83) d t a = -d t r a + l a +T a . 



For notational convenience we define the index set to be 

M = |c, i, (ij)^, bi, a | i,j = 1, 2, 3 1 . 

This set M. is just the collection of all indices in the macroscopic equations. Then 
for I 6 M. we have that each lg(t, x) are the coefficients of Z({I — P}/) with respect 
to the elements of (78); similarly each Tg(t, x) and rg(t,x) are the coefficients of 
r(/, /) and {I — P}/ respectively. Precisely, each element rg can be expressed as 

n = J2ci({l-P}f,e k ). 

k=l 

All of the constants C| above can be computed explicitly although we do not give 
their precise form herein. Each of the terms lg and Yt can be computed similarly. 

The second set of equations we consider arc the local conservation laws satisfied 
by (a^, bf ,cf). To derive these we multiply (9) by the collision invariants M in (67) 
and integrate only in the velocity variables to obtain 

d t {a f + 3c f ) + V x -b f = 
d t b f + V x (a f + 5c f ) = -V a • (v 8> vy/ji, {I - P}/) 
d t (3af + 15cf) + 5V x -bf = —V x ■ (\v\ 2 Vyfjl, {I - P}/). 

Above we have used the moment values of the normalized global Maxwcllian /i: 

(1,H) = 1, (K| a ,/*) = 1 J (M 2 ,M)=3,(|u j | 2 | Ui | 2 , M ) = l, jV<. 
(h| 4 , M )=3, (H 2 K| 2 , M ) = 5, (M 4 ,M)-15. 
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Comparing the first and third local conservation law results in 

(84) d t af = ^V a -<MV£,{I-P}/> 

(85) d t b f + V x (a f + 5c f ) = -V x -{v®v^jl,{I-P}f) 

(86) d t cf + \v x -bf = ~V x -(\v\ 2 v^JI,{I-P}f). 

3 6 

These are the local conservation laws that we will study below. For the rest of this 
section, we concentrate on a solution / to the Boltzmann equation (9). 

Lemma 10. Let f(t, x, v) be the local solution to the Boltzmann equation (9) shown 
to exist in Theorem 8 which satisfies (7). Then we have 

dx a?{t,x) = dx b (t,x) = I dx c*(t,x) = 0, 

T 3 JT 3 JT 3 

where a? , W = [&i,&2,&3], cf are defined in (68). 

The proof of this lemma follows directly from the conservation of mass, mo- 
mentum and energy (7), using the cancellation that we just used in deriving the 
conservation laws (84), (85), and (86). In the following two Lemmas, we establish 
the required estimates on the linear microscopic piece and then we estimate the 
non-linear higher order term. 

Lemma 11. For any of the microscopic terms, It, from the macroscopic equations 

l£M \a\<N 

Proof. Recall {ek]]? = i, the basis in (78). For fixed (t,x), it suffices to estimate the 
H*- 1 norm of {I ({I - P}/), e fc ). We use (77) to expand out 

(8 a l({l - P}/), e k ) = -(v ■ V X ({I - P}d a f), e k ) - (L({I - P}9 Q /), e h ). 

Now for any \a\ < N — 1 

\\(vV x ({l-P}d a f),e k )\\ 2 L2 < [ dxdv \e k (v)\ \v\ 2 |{I - P}V x d a f\ 2 

JT 3 xK 3 

<||{I-P}V^/||1^. 

Here we have used the exponential decay of e k (v). 

It remains to estimate the linear operator L. With the expression from (11) and 
Lemma 5, we have the following 

||(L({I-P}5 a /),e fc )|||, <|| \{I-P}d a f\ 1hS \M\ s J 2 Ll 

<\\{l-P}d a f\\ 2 L 2 (T 3 XR 3), (taking r, = S = 1). 

~i + 2s \ x tji 

This completes the proof of our estimates for the It. □ 

We now estimate coefficients of the higher order term /): 
Lemma 12. Let (75) be valid for some Mq > 0. Then for N > 3 we have 

E w^y? ~ E ii 9a /ii^ +2 »(TgxK3). 

eeM \ a \<N 
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Proof. As in the proof of Lemma 11, it will be sufficient to estimate the norm 
of (F(/, /), efe). We apply Lemma 5 to see that 

ll(r(/,/),e fe )|| ffN < E lli^-^/u^/UL- 

\a\<N ai<a 

We use (15), and take the supremum over the term with fewer derivatives to obtain 

< \\fh*H» E w da f\\^ +2s 

\a\<N 
|oi|<JV 

The last inequalities follow from the Sobolev embedding D □ 

We now prove the crucial positivity of L for small solution f(t, x, v) to the Boltz- 
mann equation (9). The conservation laws (7) will play an important role. 

Proof of Theorem 9. We first of all notice from (68) that 

\\Pd a f(t)\\% sn < WPaWWh + \\d a Kt)\\h x + \\d a c(t)\\k- 
Thus it will be sufficient to bound each of the terms on the right side above by 
|{I — P}d a f(t)\\fjs,-, plus the time derivative of the interaction functional, which 
is defined in (93). Indeed, our proof is devoted to establishing the following 

Mt)\\ 2 H N + \\b(t)\\i N + \\c(t)\\i N < Y ii{i-POT(*)ii!= a 

x x x £ ^ f + 2s 

\a\<N 

(87) +m E ^ a ^k +2 , + ^r- 

\a\<N 

Clearly the second term on the right above can be neglected because of 

E ii^/wiik, < E iip^/wiiI. + E ii{i-p}^/wii! 2+2 

\a\<N \a\<N \a\<N 

<{\\am H » + \\Ht)\\ H » + Mt)\\ H »} 2 + E ll{I-P}3 tt /W|||, +2s . 

\a\<N 

We have used (68). Thus (87) will imply Theorem 9 when Mq is sufficiently small. 

To prove (87), we estimate each of a, 6, and c individually with spatial derivatives 
of order < \a\ < N. Then at the end of the proof we estimate the pure I? x norm 
of a, b, and c in a uniform way. We first estimate a(t, x). Consider any |a| < N — 1. 
By taking did a of (82) and summing over i, we get 
, 3 
- Ad a a = — (V ■ 8 a b) + J2 {dtd l d a r bl - d t d a {l bl + T M }) . 

i=l 



Multiply with d a a to (88) and integrate over dx to obtain 

/ dx (V • d a b) d a a{t, x) + — [ dx did a r bi d a a{t, • 

Jf3 dt J T 3 

[ dx (V • d a b)d t d a a(t,x) - [ dx d l d a r bl d t d a a{t,x) 
+ \\d a {l bl + r bt }\\ L l\\Vd a a\\ Ll . 
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Above we implicitly sum over i = 1,2, 3. We define the interaction functional 

K (*) = / dx (V ■ d a b) d a a{t, x) + V / a 2 9 Q r fcl 5 a a(t, x). 

We also use the local conservation law (84), to see that for any 77 > 0, we have 
dx {|(V ■d a b)d t d a a{t,x)\ + \d l d a r bl d t d a a(t,x)\} 
<v\\V-d a b\\ 2 L2 +C n \\{I-P}Vd a f\\ 2 L2 . 



We combine these last few estimates with Lemma 11 and 12 to see that 
llV^oll^-^V.^Hij^C, J2 \\{I-P}d a f\\lz 2 +-£- 

\a\<N 

(89) +M l|0°7ll£» w 

|a|<iV 

This will be our main estimate for a(t, x) with derivatives. 

Next we estimate c(t,x) from (79), with \a\ < N — 1. We notice that 

||Va«c||i, < C [\\d a l c \\ 2 Ll +\\d a r c \\l i j-j t J^ dxd a r c {t,x) - V x d a c(t,x) 

+ [ dx V x ■ d a r c {t, x) d a d t c(t, x). 

We now define another interaction functional as 

= ~ I dx d a r c (t,x) ■ V x d a c(t,x). 

Next we use the conservation law (86) to obtain the following estimate 

/ dx \y x d a r c {t,x) ■d a d t c{t,x)\<i 1 \\V -d a b\\l 2 +C4{l-P}\/d a f\\l 2 , 

which holds for any rj > 0. Combining these with Lemmas 11 and 12, we see that 

dT a 

liv^ciiij-^v-^iiij^c £ \\{i-p}dv\\ 2 L2+2s + ^ 

\a\<N 

(90) +M ]T \\d a f\\ 2 L , +2s . 

\a\<N 

This will be our main estimate for c(t, x) with derivatives. 

The last term to estimate with derivatives is Vd a b. Suppose that \a\ < N — 1, 
take dj of (80) and (81) and sum on j. It was shown in a nontrivial calculation from 
[32], using the elliptic structure of these equations use several symmetries, that 

Ad a bi = -d l d l d a b i + 2d l d a k + 2d i d a T l 
-did a h ~ 1'.. + djffHij + djffTij - dtdjtiPr 
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We then multiply the whole expression by d a bi and integrate by parts to yield 

(91) \\Vd a bi\\ 2 < c J V \\d a i e \\ 2 + \\d a r e \\ 2 1 + V / dx d J d a r lJ d t d a b l 

UeM J 3 ^ Jri 

V / dx (),;Pr,,i fh,. 



d_ 

dt 



We define the last component of the interaction functional as 



With the conservation law (85), we estimate the term with a time derivative as 
J2j dx \d j d a r ij d t d a b i (t,x)\ < r, { || Vd a a\\ 2 Ll + \\Vd a c\\ 2 L ^ 

+c 1) ||{i-p}va Q /||| 2 , 



which once again holds for any 77 > 0. Combining these last few estimates with 
Lemmas 11 and 12, we obtain 

|va Q 6|| 2 L , -^{iiva-oiiij + ||va« c ||| 2 J < c, ]T ||{i-p}5«/|||, +2s + ^ 



\a\<N 



(92) 



-Mo ]T WfWl* 

|q|<At 



This is our main estimate for b(t, x) with derivatives. 

Now, with T"(t), Iu(t) and X"(t) defined just above, we define the total inter- 
action functional as 



(93) 



l(t) d = f £ {l°(t)+lZ(t)+l?(t)}. 

\a\<N-l 

Choosing for instance rj — 1/8 and collecting (89), (90), (92), we have established 

dl 

H «-x + || \irt>\\- H s-x + ||voi"c||^-i £ 2^ \W-*'W~t\\l* +2s + -fa 



\Vd a a\\ 2 rN - 1 



\Vd a bf H s- 1 + \\Vff'c\\ a H „-i< J2 \\{I-~P}d a f\ 

\a\<N 

+M £ ll^/||i, +2s . 

|a|<AT 



To finish (87), it remains to estimate the terms without derivatives. 
With the Poincare inequality and Lemma 10, a itself is bounded by 



loll < HVc 



dx 



= llVoll 



This is also bounded by the right side of (87) by the last estimate above. The 
estimates for bi(t,x) and c(t,x) without derivatives are exactly the same. This 
completes the main estimate (87) and the proof. □ 



We are now ready to prove global in time solutions to (9) exist. 
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7.3. Global Existence. With the coercivity estimate for non-linear solutions from 
Theorem 9, we prove these solutions must be global with a continuity argument. 

Proof of Theorem 1. We first fix Mq < 1 such that both Theorem 8 and 9 are valid. 
For any C > we can choose a large constant C\ > such that 

\\f(t)\\h vH » < (Ci + 1) WfWWliH- - c'T{t) < \\f(t)\\l is „. 

Notice C\ only depends upon the structure of the interaction functional and C, 
but not on f(t,x,v). We then define the equivalent instant energy functional by 

£{t)^{C 1 + l)\\f{t)\\l lH „-C'I{t). 
Then £(t) w ||/(i)||^2( K 3.jjJV(T 3 ))- Now choose Mi < ^ and consider initial data 

£(0) < Mi < M . 
From Theorem 8, we may denote T > so that 

T = sup{t > : £(t) < 2Mi} > 0. 
We now take the spatial derivatives of d a of (9) to obtain 

(94) E (Ld a f,d a f) = J2 (d a r(f,f),d a f)- 

\a\<N \a\<N 

By Lemma 2 wc have 

E (a Q r(/,/),a«/)<||/(t)|| i? H«^). 

|a|<JV 

Notice that for < t < T, by our choice of Mi, 

£{t) < 2 Ah < A1 Q . 
Thus (75) is valid. Now with Lemma 6 and then Theorem 9 we have 
E (Ld a f,d a f)>6 Q J2 ll{I-P}^|| 

\a\<N \at\<N 



2 



|a|<JV |q|<At 



{||/(*)Hi; Hj r -C'X(t)} + ~6V(t) < \\f(t)\\ LlH ,V(t). 



With 5 = min {f , ^} > and C" d = ^ > 0, we conclude that 

d_ 

dt 

We multiply (94) by C\ and add it to this differential inequality to conclude 

^ + ~8V{t)<C4f{t)\\ L 2 HS V{t), c>o. 
In the last step we have used the positivity of L > as shown in Lemma 6. Suppose 

M def . / 5 2 Mol 

iWi = mm < — — , > . 

\8cr 2 / 

We now use the definitions of Mi and T to obtain for < t < T that 
^- + 5V(t) < C4f{t)\\ LlHS V(t) < C*V2MV(t) < 5 -v(t). 
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Therefore, an integration over < t < t < T yields 

c PT 

2 Jo 

Since £{t) is continuous in r, £{T) < Mi if T < oo. This is a contradiction to the 
definition of T, thus T = oo. The time decay follows from V(t) > ||/(t)||| 2H w ■ □ 

Appendix: Carleman's representation and the dual formulation 

In this appendix we develop two Carleman [14] type representations which are 
used crucially in our main text. We consider the general expression 

C(0= / dvmv-v*\)[ da b((k,a)) H{v,v*,v',vl), 

with k = |^~^*| and the usual post-collisional velocities (t>',v*) given by (2). The 
functions b and 4> are generally given by (5) and (4). For the purposes of deriving 
the expression in Proposition 10 it suffices to suppose that both of these functions 
are smooth. The general expressions can then be deduced from these formulas by 
the usual approximation procedures. We have the following representation formula 

Proposition 10. Let H : M? x E 3 x E 3 x E 3 — > E be a smooth, rapidly decaying 
function at infinity. Then we have 

7 / / v— 2v'—v—v* \\ 

~, , f , f $(k> - vJ) ° \ \F=^T> \2v>-v-v.\ J ) 
C(«.) = 4/ dv' / dn v -iL l^^Aj \LL H 

Jr3 \v-v*\ 
Above H = H(v, v*, v' , v + u* — i/) and S£ is t/ie hyperplane 

E v v [ = f {v 1 € M 3 : («, «0 = 0} . 

TTien dvr^ denotes the Lebesgue measure on this hyperplane. 
We also derive a Carleman representation for 

C(v)= I dv* $(|«-»»|) / dab({k,a)) H(v, w„ <), 

with the same notation and the same comments as in the last case. 

Proposition 11. Lei : I 3 x K 3 x I 3 x I 3 t ie a smooth, rapidly decaying 
function at infinity. Then we have 

iff 2v — v' — v' Jf v' — u' \ \ 

V y B » * \v-v'\ \v'-v'*\ 

Above H = H(v,vl + v' — v, v', u*), and E%, is the hyperplane 

E v v , ='{<6l 3 : (t/-«X-«>=0}. 
Then dir v i denotes the Lebesgue measure on this hyperplane. 

Our expressions above may be at some degree of variance from the usual Car- 
leman representation, however they are of the same form and derived in the same 
way; a clear proof can be found in [26] . With these expressions we will derive a 
Dual Representation for the non-linear operator (10). 
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Dual Representation. Wc initially suppose that L 2 da \b((k,a))\ < oo and that 
the kernel b has mean zero, i.e., L 2 da b({k, a)) = 0. Then after the pre-post change 
of variables we can express (10) as 

(F(g,h),f)= f dv [ dv, [ da <Z>(\v - v*\)b ((k,a)) g*h (M'J' - M,f) 
JR3 JM. 3 Js 2 



dv / dv, / da^{\v- Vlt \)b((k,<T))g t hM'J'. 
This follows from the vanishing of J g2 b((k, a))da. With Proposition 10, this is 

= i[dvj dv' I dn v *(|« - ^|) & (\ 1 7^ T ' ^ ^H g .hM'J'. 



As usual above d7r„ is Lebesguc measure on the two-dimensional plane E% passing 
through v 1 with normal v' — u*, and of course v is the variable of integration. In 
the above formulas, we take Ml = M(v + u* — v'). From the identity 



v — v* 2v' — v—v* \ \v' — v*| 2 — \v — v'\ 



l\2 



\v — v*\ ' \2v' — v — v*\ I \v — v'\ 2 + \v' — v*| 2 ' 
we observe that 

[ a hi I v ~ v * 2v> ~ v ~ v * \\ K ~ v *\ 2 

J k: dn ' v [\\v- v.\' \2v> -v-vA/J \v- 

Jo \\v '-v*\ 2 +r V (r 2 + \v '-v*\ 2 ) 

by a change of variables since J_ x di 6(i) = and 



d_ 

dr 



\v' — vj 2 — r 2 



-4r\v' - v*\ 2 
(r 2 + \v' - v*| 2 ) 2 ' 



\v' — v, \ 2 + r 2 
In particular, this implies 

I dn ZJW - / /v-v. 2v> - v - v \\ kV-^ /Jf = Q 

We subtract this expression from the Carlcman representation just written for 
(T(g, h), /), to see that (T(g, h),f) must also equal 



ev'\v' — v*\\v — u*| \\ \ v ~ v *\ \^ v ' — v — v*\ 
x ($(v- vJhMi - - «.) ^ ~ h'M?\ . 

V \v-vA 3 J 



This will be called the "dual representation." 

The claim is now that this representation holds even when the mean value of the 
singular kernel b((k, a)) from (4) is not zero. To see this claim, suppose b integrable 
but without mean zero. Then define 



r-l 

-1- 



b t {t) = b(t)-€- 1 l [1 ^ 1] {t) I b{t)dt. 
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As a function on § 2 , b e will clearly have a vanishing integral. However, given 
arbitrary /, g and h which are Schwartz functions, it is not hard to see that 

(T(g,h),f)-(T e (g,h)J)^0, e -> 0. 

Above r e is the non-linear term (10) formed with b € (t) in place of b(t). This 
convergence holds because cancellation guarantees that the integrand vanishes on 
the set defined by (k,a) = 1. Moreover, an additional cutoff argument shows that 
the equality also holds provided that b(t) satisfies (4); the higher-order cancellation 
is preserved because ^ v _"*j possesses radial symmetry in v — v' . 
The "dual representation" deserves its name because if one defines 

T g f(v) = / dv* f daBg* {M'J' - M*f) 

W)^/ dv J -. *± b (/f^ "-"* )) g . 

JK3 J e£ \v — V *\\ v v *\ \\\ v ~ v *\ \2v'-v-v*\/J 

x - v*)Mlh - - «.) Y ~ V *l Mth 1 ) , 

then 

(95) (T(g,h),f) = {T g f,h) = (f,T;h). 

Note that the last inner product above represents an integration over dv' whereas 
the first two inner products above represent integrations over dv. 

The advantage of this representation is that T g f and T*h both depend on g in 
a fairly elementary way. This will allow, for example, the trilinear form (T(g, h), /) 
to be understood as a superposition of bilinear forms in h and /. 
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